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The  goal  of  this  research  program  was  a  quantitative  undent  ending  of  turbulence -chemistry  in¬ 
teractions  pertinent  to  future  aeropropulsion  combustors.  For  example,  (t)  flameout  and  relight  in  tur¬ 
bine  combustors  are  related  to  interactions  of  turbulence  with  chain-branching  reactions,  (ii)  hydrogen 
burnout  in  supersonic  combustors  is  related  to  interactions  with  recombination  reactions,  and  (iii)  em¬ 
issions  of  NO,,  CO,  smoke  and  other  observables  are  related  to  noocquiUbrium  in  the  populations  of 
intermediate  species  such  as  oxyhydrogen  radicals  and  C,Hy. 

A  bluff-body  stabilized  turbulent  diffusion  flame,  time-  and  space-resolved  laser  Raman  measure¬ 
ments  of  major  species,  and  a  nonequilibrium  computational  fluid  mechanics  code  were  applied  to  the 
problem.  The  principal  conclusions  art 

t  Turbulent  jet  flames  are  being  abandoned  in  the  search  for  more  intensely  turbulent  flames. 

1  An  axarymmetric  bluff-body  stabflbied  turbulent  diffusion  flue  burner  is  a  reasonable  choice  for 
combustion  research  at  hi^i  Reynolds  numbers,  approaching  btowott. 
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3.  Raman  scattering  for  measurements  of  major  spedes  and  temperature  can  be  extended  into  the 
sooting/chemiluminescent  environment  of  methane  flames.  Space-  and  time-resolved  Raman 
scattering  measurements  have  been  made  in  bluff-body  stabilized  CO/H2/N2  and  CH*  flames  at 
conditions  approaching  blowoff. 

4.  A  thermochemical  submodel  based  on  partial  equilibrium  in  the  oxyhydrogen  radical  pool  was 
developed  for  the  27.5%  CO/  323%  H2/40 2%  N2-air  system.  The  chemistry  can  be  described 
in  terms  of  two  scalars.  The  elliptic  form  of  the  time-averaged  Navier  Stokes  equations  with  k-e 
closure  was  solved  using  an  iterative  finite-volume/pressure-correction  algorithm.  Mean  pro¬ 
perties  such  as  density  were  obtained  at  each  iteration  by  convolution  with  the  joint  pdf  over  the 
two  thermochemical  scalars.  The  pdf  itself  was  obtained  by  the  moment-equation/assumed- 
shape  method. 

5.  Calculations  compare  favorably  with  the  Raman  data. 

6.  Work  on  pdf  transport/Monte-Carlo  methods  in  recirculation-stabilized  flames  has  begun. 

7.  The  range  of  turbulence-chemistry  interactions  in  combustion  has  been  quantified,  in  an  attempt 
to  gauge  the  universality  of  turbulence-chemistry  models.  While  some  combustion  reactions  fall 
in  the  flamelet  regime  and  some  in  the  distributed  regime,  many  important  reactions  fall  in 
between.  Turbulence-chemistry  interactions  occur  in  a  greater  variety  of  modes  than  accounted 
for  by  current  combustion  theory.  The  simultaneous  prediction  of  kinetically  sensitive  phenome¬ 
na  such  as  pollutants,  flame  stability  and  combustion  efficiency,  will  require  a  new  approach. 
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Section  1 
SUMMARY 

This  report  describes  work  performed  in  the  period  May  1988  to  June  1991  for  the  U.S.  Air  Force 
Office  of  Scientific  Research  under  Contract  F49620-88-C-0066.  Dr.  Julian  Tishkoff  was  the  AFOSR 
Project  Monitor.  The  work  was  performed  and  the  report  has  been  prepared  by  the  Engineering  Sys¬ 
tems  Laboratory  at  the  General  Electric  Research  and  Development  Center,  Schenectady,  New  York. 
This  section  summarizes  the  technical  problem  which  was  considered  and  outlines  the  results  ob¬ 
tained  in  the  course  of  the  program. 

Aircraft  and  industrial  gas-turbine  engines  generally  employ  non-premixed  turbulent  combustion. 
The  goal  of  this  three-year  research  program  was  a  quantitative  understanding  of  turbulence- 
chemistry  interactions  in  such  flames.  Although  the  fast  chemistry  or  “mixed-is-burned”  paradigm  is 
useful  for  understanding  traditional  design  issues,  present-day  demands  on  combustion  equipment 
and  thus  on  computational  models  are  increasing  the  need  for  a  more  advanced  understanding  of 
turbulence-chemistry  interactions.  For  example,  (i)  flameout  and  relight  in  turbine  combustors  are 
related  to  interactions  of  turbulence  with  chain-branching  reactions,  (ii)  hydrogen  burnout  in  super¬ 
sonic  combustors  is  related  to  interactions  with  recombination  reactions,  and  (iii)  emissions  of  N00 
CO,  smoke  and  other  observables  are  related  to  nonequilibrium  in  the  populations  of  intermediate 
species  such  as  oxyhydrogen  radicals  and  QHr 

Statement  of  Work 

a.  Identify  conditions  necessary  to  observe  extinction  in  pilot-stabilized  hydrocarbon  jet  flames. 

b.  Modify  the  Raman  system  to  obtain  data  in  the  blue  region  of  hydrocarbon  jet  flames;  measure 
major  species  concentrations,  temperature,  and  density  under  local  extinction  conditions;  meas¬ 
ure  inlet  conditions  for  modeling  purposes. 

c.  Extend  the  chemical  scheme  in  the  model  to  hydrocarbons,  using  multistep/partial-equilibrium 
schemes  and  apply  joint  pdf/Monte  Carlo  and  conventional  flow  models  to  hydrocarbon  jet 
flames. 

d.  Compare  data  and  calculations;  evaluate  the  model  and  mechanisms  of  extinction,  hypothesize 
new  mechanisms  if  necessary. 

e.  Use  a  combustor  with  bluff-body  stabilization  of  carbon  monoxide/hydrogen/nitrogen  flames  at 
high  Reynolds  numbers  and  establish  local  extinction  conditions. 

f.  Measure  major  species  concentrations,  temperature,  and  density  in  bluff-body  stabilized  carbon 
monoxide/hydrogen/nitrogen  flames  under  conditions  of  significant  local  extinction;  if 
significant  levels  are  not  found,  use  hydrocarbon  fuels;  measure  inlet  conditions  for  modeling 
purposes. 

g.  Model  the  bluff-body-stabilized  flames  using  state-of-art  numerical  techniques,  algebraic  or 
differential  Reynolds  stress  models  for  turbulence,  and  turbulence-chemistry  interaction  models 
suggested  by  the  work  in  jet  flames,  compare  with  data  and  draw  conclusions  regarding  extinc¬ 
tion  mechanisms. 

h.  In  flames(s)  that  show  the  most  significant  local  extinction,  measure  velocity  and  turbulence  in¬ 
tensity  (LDV)  and  minor  species  concentrations  (saturated  laser-induced  fluorescence  of  OH), 
establish  the  topology  of  the  flame  through  planar  OH  imaging  and  Schlieren  photography,  esti¬ 
mate  scalar  dissipation  rates  by  Rayleigh  (or  Raman,  if  possible)  imaging. 
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i.  Investigate  the  universality  of  the  turbulence-chemistry  interaction  model  by  calculating  a  range 
of  measured  nonequilibrium  phenomena  from  low  levels  of  superequilibrium  radicals  to  local 
extinction. 

Status 

A  bluff-body  stabilized  turbulent  diffusion  flame,  time-  and  space-resolved  laser  Raman  measure¬ 
ments  of  major  species,  and  a  nonequilibrium  computational  fluid  mechanics  code  were  applied  to  the 
problem.  The  principal  conclusions  are 

1.  Turbulent  jet  “diffusion”  flames  are  being  abandoned  in  the  search  for  more  intensely  turbulent 
flames.  Jet  flames,  however,  have  had  a  highly  successful  record  over  the  last  twenty  years  and 
have  provided  most  of  the  current  understanding  of  turbulence-chemistry  interactions  in  non- 
premixed  flames. 

2.  An  axisymmetric  bluff-body  stabilized  turbulent  diffusion  flame  burner  is  a  reasonable  choice 
for  combustion  research  at  high  Reynolds  numbers,  approaching  blowoff. 

3.  Raman  scattering  for  measurements  of  major  species  and  temperature  was  extended  into  the 
sooting/chemiluminescent  environment  of  methane  flames.  (This  requires  substantial  time  and 
effort  and  precluded  LDV  and  LIF  (OH)  measurements.)  Space-  and  time-resolved  Raman 
scattering  measurements  have  been  made  in  bluff-body  stabilized  CO/H2/N2  and  CH4  flames 
at  conditions  approaching  blowoff. 

4.  A  thermochemical  submodel  based  on  partial  equilibrium  in  the  oxyhydrogen  radical  pool  was 
developed  for  the  27.5%  CO/  32.3%  H2 /  40.2%  N2-air  system.  The  chemistry  can  be 
described  in  terms  of  two  scalars. 

5.  The  elliptic  form  of  the  time-averaged  Navier  Stokes  equations  with  k-e  closure  was  solved  us¬ 
ing  an  iterative  finite-volume/pressure-correction  algorithm.  Mean  properties  such  as  density 
were  obtained  at  each  iteration  by  convolution  with  the  joint  pdf  over  the  two  thermochemical 
scalars.  The  pdf  itself  was  obtained  by  the  moment-equation/assumed-shape  method. 

6.  Calculations  have  been  compared  with  the  Raman  data. 

7.  Work  on  pdf  transport/Monte-Carlo  methods  in  recirculation-stabilized  flames  has  begun. 

8.  The  range  of  turbulence-chemistry  interactions  in  combustion  has  been  quantified  in  an  attempt 
to  gauge  the  universality  of  turbulence-chemistry  models.  While  some  combustion  reactions  fall 
in  the  flamelet  regime  and  some  in  the  distributed  regime,  many  important  reactions  fall  in 
between. 

This  study  has  shown  that  turbulence-chemistry  interactions  occur  in  a  greater  variety  of  modes 
than  accounted  for  by  current  combustion  theory.  The  simultaneous  prediction  of  kinetically  sensitive 
phenomena  such  as  pollutants,  flame  stability  and  combustion  efficiency  will  require  a  new  approach. 


2 


Section  2 

TECHNICAL  DISCUSSION 

This  study  consisted  of  experiments  on  and  modeling  of  turbulence-chemistry  interactions  in  tur¬ 
bulent  gaseous  flames  under  conditions  approaching  blowoff.  It  builds  on  prior  work,  first  on  unas¬ 
sisted  and  more  recently  on  pilot-stabilized  CO/H2/N2  jet  flames  in  coflowing  air.  To  obtain  higher 
turbulent  strain  rates  than  attainable  even  with  pilot-flame  stabilization,  an  axisymmetric  bluff-body 
stabilized  burner  was  used.  First,  however,  the  merits  of  jet  flames  are  reviewed. 

The  relatively  simple  fuel-jet-in-coflowing-air  configuration  has  been  favored  for  twenty  years  be¬ 
cause  the  mean  profiles  can  be  predicted  by  gradient  diffisiop.  models  [Starner  and  Bilger  (1980, 
1981);  Lockwood  and  Stolakis  (1983)],  the  flow  can  be  described  by  the  computationally  straightfor¬ 
ward  shear  layer  form  of  convection-diffusion-source  equations,  and  the  flame  permits  laser-based 
spectroscopic  measurements  which  have  become  critical  to  progress  in  understanding.  Penner  et  al. 
(1984)  review  a  wide  variety  of  these  laser-based  methods.  Jet  flames  can  be  fully  turbulent  and  man¬ 
ifest  a  wide  range  of  turbulence-chemistry  interactions.  Because  of  these  advantages,  studies  in  jet 
flames  have  contributed  much  to  turbulent  combustion. 

As  the  Reynolds  number  of  the  jet  is  increased,  the  flame  blows  off  the  lip  at  the  exit  of  the  nozzle. 
Thus  strong  turbulence-chemistry  interactions  leading  10  turbulent  extinction  cannot  be  studied. 

In  the  preceding  study,  we  investigated  a  non-premixed  turbulent  jet  flame  composed  volumetri- 
cally  of  a  40%  CO/10%  H2/50%  N2  mixture  in  co-flowing  air  [Correa  and  Gulati  (1989)].  This  mix¬ 
ture  of  CO  and  H2  was  used  as  the  fuel  because  the  chemistry  is  known  and  can  be  simplified  reason¬ 
ably  well  and  the  extinction  limits  approach  those  of  hydrocarbons,  which  are  the  fuels  of  practical 
relevance.  The  flame  was  stabilized  by  a  co-annular  premixed  pilot  flame.  Pointwise  measurements 
included  (i)  temperature,  mole  fractions  of  major  species,  and  mixture  fraction  £,  by  spontaneous  Ra¬ 
man  spectroscopy  and  (ii)  velocity,  by  laser  yelocimetiy.  The  data,  e.g.,  temperature-mixture  fraction 
scattergrams,  showed  significant  temperature  decrements  due  to  finite-rate  chemistry  but  no  evidence 
of  localized  extinction. 

A  computational  model  for  a  jet  flame  under  such  conditions  was  also  presented  in  Correa  and 
Gulati  (1989).  Combustion  chemistry  was  represented  by  two-body  shuffle  reactions  and  three-body 
recombination  reactions.  The  scalar  dissipation  rate  field  was  examined  for  a  critical  value  below 
which  the  two-body  reactions  were  assumed  to  be  in  partial  equilibrium  and  above  which  they  were 
assumed  to  be  frozen  and  the  gas  therefore  unburned.  The  kinetics  of  the  recombination  reactions 
were  activated  for  the  former  fraction  of  the  gas.  This  approach  was  implemented  in  a  shear-layer 
finite-volume  averaged  Navier-Stokes  model  with  k-e/assumed  shape  probability  density  function 
(pdf)  submodels  for  turbulence.  It  was  concluded  that  local  strain-induced  extinction  is  not  seen  in 
the  above  CO/H2/N2-air  flame,  despite  its  relatively  low  critical  strain  of  950  s'1  in  a  laminar 
counterflow  diffusion  flame  (“Tsuji”)  burner.  Localized  extinction  is  absent  because  the  reaction 
zones,  centered  on  ^cs0.43,  are  too  broad  to  interact  with  the  straining  scales  [Peters  ( 1988),  Bilger 
( 1989)].  H2  and  CH4  flames  have  much  narrower  reaction  zones,  and  CH4  flames  do  indeed  show  lo¬ 
calized  extinction.  Alkyl  consumption  of  radicals  may,  however,  provide  a  chemical  (nonaerodynam- 
ic)  explanation  for  extinction  in  the  methane  case  [Bilger  (1989)]. 

The  jet-flame  work  has  shown  that  it  is  overly  simplistic  to  consider  turbulent  flames  as  ensembles 
of  laminar  flames  subjected  to  strain.  The  reason  is  that  different  kinetic  steps  in  a  combustion 
mechanism  can  have  widely  disparate  time-scales  (e.g.,  chain-branching  reactions  are  much  faster 
than  recombination  reactions).  The  fastest  reactions  can  interact  with  the  local  turbulence  in  a 
“flamelet”  manner  while,  for  the  slow  reactions,  the  Damkohler  numbers  are  relatively  small  and  the 
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flame  is  more  akin  to  a  well-distributed  zone  of  reaction.  Section  2.7  presents  a  more  detailed  discus¬ 
sion  of  this  point. 

2.1  Choice  of  Axisymmetric  Bluff-Body  Stabilized  Burner  over  Jet  Flame 

To  increase  the  interaction  between  turbulence  and  chemistry,  and  the  probability  of  localized  ex¬ 
tinction,  the  jet  velocity  in  the  pilot-stabilized  burner  was  increased;  however,  only  a  small  increase 
was  possible  before  the  flame  blew  off  the  lip.  Moreover,  the  coannular  premixed  pilot  flame  used  to 
stabilize  the  (non-premixed)  jet  flame  complicates  the  near-field  turbulence,  without  greatly  increas¬ 
ing  permissible  dissipation  rates  prior  to  blowoff.  To  remedy  this,  we  switched  to  an  axisymmetric 
bluff-body  stabilized  non-premixed  flame,  with  a  view  to  obtaining  Raman  data  or  turbulence- 
chemistry  interactions  at  very  high  dissipation  rates.  While  axisymmetric  bluff-body  stabilized  burners 
can  be  operated  in  a  regime  dominated  by  unsteady  vortex  shedding  [Correa  and  Pitz  (1982),  Ro- 
quemore  et  al.  (1982)],  here  we  wish  to  avoid  unsteadiness  in  the  mean  flow. 

Axisymmetric  bluff-body  stabilized  flames  have  two  advantages  over  jet  flames,  but  one  significant 
drawback.  The  first  advantage  is  that  the  cylindrical  shear  layer  behind  the  bluff  body  provides  a  re¬ 
gion  of  large  dissipation  rates,  located  far  enough  downstream  of  the  bluff  body  that  perturbations  at 
the  inlet  should  not  affect  the  dissipation  field  dramatically.  [In  jet  flames  the  maximum  dissipation 
rates  are  at  the  lip  of  the  jet,  a  region  which  is  very  sensitive  to  the  inlet  conditions.]  The  second  ad¬ 
vantage  is  that  there  is  no  need  for  pilot  stabilization  and  thus  there  is  no  interference  from  the 
premixed  pilot  as  might  occur  in  jet  flames.  On  the  other  hand,  the  major  drawback  of  bluff-body 
flows  is  that  they  cannot  be  modeled  as  accurately  as  shear  flows  and  jets  because  of  the  elliptic  na¬ 
ture  of  the  flowfield  and  the  lack  of  sufficiently  reliable  turbulence  models  for  recirculating  flow. 
These  problems  historically  overshadowed  the  modeling  of  many  turbulence-chemistry  effects,  mak¬ 
ing  evaluation  of  the  models  difficult.  Large  supercomputers  and  improved  modeling  techniques, 
however,  have  made  it  worthwhile  to  model  the  axisymmetric  bluff-body  experiments  now  being  per¬ 
formed. 

Other  researchers  are  also  being  driven  away  from  jet  flames  in  the  search  for  higher  turbulence 
levels.  Notice  the  increase  in  complexity  in  the  following  sequence  of  experiments: 

1.  Simple  jet  flames,  e.g.,  those  reviewed  by  Strahle  and  Lekoudis  (1985). 

2.  Pilot-stabilized  jet  flames  (pilot  disturbances  in  enthalpy  and  other  quantities  are  generally  as¬ 
sumed  away  in  calculations). 

3.  The  axisymmetric  “Reverse  Flow  Reactor,”  employed  at  the  University  of  Sydney/SANDIA 
(Mansour  et  al.  1989),  which  uses  flow  recirculation  to  stabilize  the  flame  above  jet-blowoff  con¬ 
ditions. 

4.  Axisymmetric  bluff-body  stabilized  flames,  such  as  one  being  studied  using  advanced  laser  diag¬ 
nostics  including  Raman  imaging  by  an  ALTEX/SANDIA  team  [Schefer  et  al.  (1989)].  An  ear¬ 
lier  propane-fueled  bluff-body  stabilized  flame  at  the  U.S.  Air  Force  Aeropropulsion  Laborato¬ 
ry  [Roquemore  et  al.  (1982)]  was  analyzed  computationally  by  Correa  (1984),  using  an  equilibri¬ 
um  chemistry/assumed  shape  pdf  approach.  This  study  revealed  many  of  the  shortcomings  of 
such  models.  Later  work  by  Switzer  et  al.  (1985)  applied  CARS  to  this  bluff-body  stabilized 
flame. 

5.  The  GE  axisymmetric  bluff-body  stabilized  burner  developed  and  used  in  the  present  study. 

It  is  clear  that  simple  “parabolic”  flows  are  being  replaced  in  the  drive  toward  strong  turbulence- 
chemistry  interactions.  It  will  be  necessary  to  accurately  model  these  more  complex  flows  if  the  link 
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between  theory  and  experiment  is  to  continue  as  productively  as  before. 

22  Setup  of  Experiments 

Since  the  turbulence  in  a  pilot-stabilized  non-premixed  turbulent  CO/H2/N2  jet  flame  was  not  in¬ 
tense  enough  to  cause  local  extinction,  axisymmetric  bluff-body  stabilization  was  resorted  to. 

22.1  Bluff-Body  Stabilized  Burner 

Typical  conditions  and  a  schematic  of  the  burner  in  the  tunnel  are  shown  in  Figure  1.  The  subson¬ 
ic,  low  turbulence  wind  tunnel  has  been  used  extensively  for  experimental  studies  of  laminar,  transi¬ 
tional,  and  turbulent  diffusion  flames  in  the  past  [Lapp  et  al.  (1983)].  The  axisymmetric  bluff-body  has 
an  outer  diameter  of  38.1  mm  with  an  axial  jet  of  diameter  3.18  mm  located  in  the  center.  The  back 
surface  of  the  bluff  body  is  coated  with  a  thermal  barrier  material  to  reduce  heat  loss.  The  flame  is 
stabilized  by  the  recirculation  zone  provided  by  the  bluff  body.  The  ratio  of  the  bluff-body  diameter  to 
the  central  jet  diameter  is  12:1.  The  tunnel  cross  section  is  6  in.  x  6  in.  so,  and  the  bluff-body  blockage 
is  -  5%.  Comprehensive  observations  and  photographs  of  H2,  CO/H2/N2,  CH4,  and  CH4/N2 
flames  at  various  air  and  fuel  jet  velocities  were  obtained  earlier  and  have  been  used  to  select  promis¬ 
ing  cases  for  the  Raman-scattering  measurements  repot '  below. 

222  Raman  Setup 

The  Raman  diagnostic  system  used  for  space  (almost)  and  time-resolved  measurement  of  tem¬ 
perature,  mixture  fraction  and  major  species  is  similar  to  that  used  earlier  for  H2  and  H2/CO/N2 
flames  and  was  described  in  detail  by  Lapp  et  al.  (1983).  For  Raman  measurements  in  methane 
flames,  the  system  was  extensively  modified  to  account  for  the  significant  levels  of  laser-induced 
fluorescence/incandescence  encountered  in  the  rich  sooty  zone  of  flame,  as  described  below. 

(a)  Raman  Measurements  in  Bluff-Body  Stabilized  CO/H2/N2  Flame 

The  Raman  system  is  based  on  a  flashlamp  pumped  dye  laser  which  provides  pulses  of  ~  1J  in  2-4 
/is  within  a  0.2  nm  bandpass  at  488.0  nm  at  a  repetition  rate  of  10  pps.  The  light  scattered  at  right  an¬ 
gles  is  collected  by  two  lenses,  separated  in  frequency  by  a  3/4  m  spectrometer  and  is  detected  by 
eight  photomultiplier  tubes.  The  photomultiplier  tubes  detect  anti-Stokes  vibrational  Raman  scatter¬ 
ing  from  N2,  Stokes  vibrational  Raman  scattering  from  N2,  02,  H2,  H20,  CO  and  C02,  and  Rayleigh 
scattering.  The  temporal  resolution  (2  /us)  of  the  technique  is  limited  by  the  laser  pulse  length,  the 
spatial  resolution  (0.2  x  0.2  x  0.6  mm)by  the  spectrometer  entrance  slit  and  the  collection  optics,  while 
the  data  aquisition  rate  (5  Hz)  is  limited  by  the  laser  repetition  rate.  The  flame  luminescence  was 
found  to  be  very  broad-banded  throughout  the  visible  region.  A  polarization  filter  was  used  in  the 
collection  optics  system  to  reduce  flame  luminescence.  The  polarization  vector  was  aligned  to  pass 
the  horizontally  polarized  Raman-  and  Rayleigh-scattered  light.  Typically  200  measurements  were 
made  at  each  flame  location.  2000  shots  were  also  recorded  at  some  locations  for  statistical  purposes. 

The  joint  Rayleigh- Raman  data  is  obtained  as  follows.  The  photomultiplier  tubes  are  gated  elec¬ 
tronically  before  and  after  the  firing  of  the  laser.  The  differential  voltages,  which  correspond  to  the 
Raman  signal  minus  the  flame  background  at  t  =  0  are  amplified  and  filtered  before  digitization  using 
a  12-bit  A/D  convertor.  The  signals  are  then  corrected  for  electrical  and  other  background  errors 
corresponding  to  each  channel  and  are  normalized  by  the  laser  energy  for  every  shot.  The  subsequent 
signals  are  corrected  for  relative  sensitivities  of  the  photomultiplier  tubes  to  obtain  the  final  values 
corresponding  to  the  vibrational  intensities  of  Raman  scattering  (which  is  linearly  proportional  to  the 
number  density)  in  the  bandpass  of  the  exit  slits  of  the  spectrometer.  The  system  is  calibrated  exten- 
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sively  using  100%  pure  gases  and  well-characterized  premixed  porous  plug  burners  before  and  after 
each  set  of  data  in  the  turbulent  flame  is  obtained.  The  data  is  collected  using  data  translation 
DT2821  board  mounted  in  a  DELL333  system  using  GLOBAL  lab  software-.  The  data  is  processed 
on  the  IBM  compatible  DELL  system.  Figure  2  shows  the  Raman  system  consisting  of  the  dye  laser, 
data  acquisition  electronics,  and  the  flow  tunnel.  Figure  3  shows  the  spectrometer  and  the  PMT 
housing  box  from  the  other  side  of  the  setup. 

A  photograph  of  the  bluff-body  stabilized  CO/H2/N2  flame  used  in  this  study  is  shown  in  Figure 
4.  The  flame  is  bluish-white  because  of  chemiluminescence,  but  no  yellow  sooty  region  is  visible. 

The  instantaneous  temperature  at  every  shot  was  determined  using  three  independent  methods. 
The  Stokes-anti-Stokes  (SAS)  ratio  method  as  described  in  Drake  et  al.  (1982)  is  based  on  the  nitro¬ 
gen  element  and  yields  the  temperature  directly.  The  second  method  is  based  on  an  iterative  scheme 
in  which  an  initial  temperature  is  guessed,  based  on  which  the  mole  fraction  of  all  major  species  are 
calculated  using  their  measured  vibrational  intensities.  The  mole  fractions  are  then  corrected  using 
high-temperature  correction  factors  to  account  for  changes  in  the  fraction  of  the  Raman  band  falling 
in  the  exit  slits  provided  for  the  respective  photomultiplier  tubes.  The  iteration  process  is  repeated  till 
the  sum  of  the  mole  fractions  is  unity.  This  iteration  procedure  converges  in  two  to  three  iterations 
since  the  correction  factors  are  relatively  weak  functions  of  temperature.  The  mixture  fraction  is  cal¬ 
culated  using  the  definition  suggested  by  Bilger  (1989)  which  takes  into  account  the  balance  of  three 
atoms  i.e.,  C,H,  and  O,  and  is  defined  such  that  the  conserved  scalar  preserves  its  stoichiometric  value 
even  in  the  presence  of  differential  diffusion. 


2  Zc/Wc  +  rrZH/w„  +  (Z0,„  -  z0)/w0 


2  ZQf/Wc  +  i  Z„j/WH  +  Zo,JW0 

Other  definitions  of  mixture  fraction  were  also  used. 

In  this  flame  it  was  found  that  the  chemiluminescence  caused  by  the  CO  +  O  -*•  C02  reaction  was 
significant  and  extra  precautions  had  to  be  taken  to  improve  the  signal-to-noise  ratio,  such  as  inserting 
a  polarization  filter  in  the  collection  optics,  reducing  the  entrance  slit  size,  and  so  on.  In  addition, 
since  the  N2  channel  is  affected  most  by  the  luminosity  of  the  flame,  the  SAS  measurements  were  sus¬ 
ceptible  to  this  error.  Hence,  Rayleigh  scattering  was  used  to  obtain  an  independent  measurement  of 
temperature.  Since  the  flame  is  of  the  diffusion  or  non-premixed  type,  additional  information  about 
the  Rayleigh  cross  section  of  the  species  involved  is  needed  for  Rayleigh  scattering  to  be  applicable. 

In  this  flame  the  Raman  scattering  signals  from  the  major  species  were  used  to  obtain  information 
about  changes  in  the  Rayleigh  cross  section  of  the  mixture,  and  thus  we  were  able  to  measure  tem¬ 
perature  in  an  iterative  manner. 

Results  are  discussed  in  subsection  2.3. 

(b)  Raman  Measurements  in  a  Bluff-Body  Stabilized  Pure  Methane  Flame 

Figure  5  shows  that  the  methane  flame  has  a  sooty  yellow  region  in  the  recirculation  zone  of  the 
bluff  body.  The  flame  then  narrows  significantly  and  appears  to  be  intermittently  connected  in  the 
neck  region,  following  which  the  flame  is  fairly  blue.  Initial  measurements  with  the  Raman  system 
showed  that  there  was  significant  laser-induced  fluorescence  interference  throughout  the  flame,  as 
also  reported  in  Masri  et  al.  (1987, 1988, 1989),  and  Dibble  et  al.  (1987).  The  fluorescence  is  believed 
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Figure  2.  Photograph  of  the  Raman  system  consisting  of  the  dye  laser,  data 
acquisition  electronics,  and  the  flow  tunnel 


I 

f  Figure  3.  Photograph  shewing  the  spectrometer  snd  the  PMT  housing  box. 

The  bluff-body  has  a  blockage  of  5%. 
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to  be  from  soot  precursors  and  other  higher  hydrocarbons.  Profiles  suggested  that  the  laser-induced 
fluorescence  was  fairly  broadband  and  contaminated  all  Raman  channels,  but  surprisingly  was  negligi¬ 
ble  in  the  Rayleigh  channel.  An  additional  problem  encountered  in  the  application  of  Raman  to  this 
flame  is  the  crosstalk  between  CH  and  otu0*-  Raman  channels,  primarily  02.  For  the  488.0  nm  excita¬ 
tion  used  here,  the  Raman  interference  ir.  i-  major  species  was  insignificant. 

To  account  for  these  two  additional ' .  .  s  f  contamination  in  the  signals,  the  system  as  well  as 

calibration  procedures  were  signific  *■  V  ■nj.';  ,..J  ‘a  the  following  manner.  Since  the  laser-induced 
incandescence/fluorescence  interference  :  „  ..ty  broadband,  additional  photomultiplier  tubes  were 
installed  in  two  Raman-free  regions  of  the  ; .  .ectra  to  monitor  the  laser-induced  fluorescence  on  a 
shot-to-shot  basis.  Figure  6  shows  the  locat  an  of  these  two  additional  PMTs  in  the  box  relative  to  the 
existing  photomultiplier  tube  locations.  Th.  ,  r'MTs,  termed  FI  and  F2  respectively,  were  located  at 
540  nm  (between  02  and  C02)  and  at  590  >  a  (between  CH4  and  H20)  respectively.  These  two  sig¬ 
nals  were  found  to  correlate  very  we’i  with  .  acn  other  and  with  all  other  Raman  signals,  so  that  in  the 
calibration  procedures  used  eventually,  the  use  of  FI  was  found  to  be  sufficient  to  allow  corrections  in 
all  other  channels.  Since  space  was  severely  limited  in  the  PMT  housing  box,  the  new  PMTs  were 
mounted  remotely  with  a  fibre  optiCo  bundle  connecting  the  slit  exits  to  the  cathodes  of  the  PMTs. 

Tne  calibration  procedure  used  to  correct  for  laser  induced  fluorescence  is  similar  to  that  used  by 
Dibble  et  al  (1987, 1990).  To  facilitate  daily  calibrations,  a  38.0  mm  dia.  honeycomb  burner  was  built 
to  provide  laminar  diffusion  flames  of  CH4  and  CO  stabilized  at  the  burner  end.  The  burner  had 
flame  arrestors  and  flow  straightners  in  it  and  resulted  in  a  smooth  thin  laminar  flame  which  was  visi¬ 
bly  sooty  and  yellow  at  the  downstream  end,  but  contained  enough  soot  precursors  and  laser-induced 
Ouorescence  at  the  upstream  end  to  allow  calibrations.  The  amount  of  soot  and  the  color  of  the  flame 
could  be  easily  controlled  by  varying  the  ratio  of  CH4  to  CO  in  the  laminar  flame.  Figures  7  and  8 
present  profiles  of  Raman  data  for  major  species  obtained  10  mm  downstream  of  the  burner  exit  in  a 
laminar  CH4/CO  (70%)  fame.  Measured  laser-induced  fluorescence  signals  monitored  /ia  channel 
FI  are  also  plotted.  On  the  r.ch  side  of  the  flame  contamination  of  Raman  signals  is  observed  in  all 
major  species  as  evident  by  the  “bumps”  in  the  raw  data  profiles.  The  shape  of  the  bump 
corresponds  to  the  shape  of  the  FI  signal  for  all  channels. 

Additional  tests  showed  that  the  magnitude  of  the  bump  also  increased  linearly  with  the  amount  of 
soot  in  the  flame,  and  with  increased  FI.  For  major  species  such  as  N2,  CH4,  and  CO  the  contamina¬ 
tion  due  to  laser-induced  fluorescence  can  be  readily  accounted  for  because  they  are  expected  to  de¬ 
cay  monotonically  through  the  flamefront.  The  situation  for  O*  and  product  species  such  as  C02  and 
H20,  is  more  complicated.  02  is  also  contaminated  by  cross-talk  from  the  CH4  channel.  The  pro¬ 
duct  species  C02  and  H20  peak  somev/here  in  the  flamefront,  where  FI  also  peaks. 

The  calibration  factors  to  correct  for  fluorescence  in  all  channels  were  obtained  by  iteration.  First 
the  raw  data  were  used  (without  an)  correction  factors)  to  calculate  temperature  (using  the  sum  of 
mole  fraction  method)  and  mole  fraction  of  major  species.  The  calculated  results  were  plotted  in  mix¬ 
ture  fraction  space,  and  compared  with  equilibrium  and  counterflow  laminar  diffusion  flame  calcula¬ 
tions  for  methane  in  air  under  assumed  values  of  stretch  [Chen  et  al.  (1989)].  Figures  9  to  11  show 
typical  results.  The  effect  of  contamination  caused  by  fluorescence  is  clearly  evident  from  the  kinks  in 
the  measured  profiles,  particularly  in  the  fuel-rich  regions  of  the  flame.  The  next  step  in  the  process 
is  to  estimate  correction  factors  for  each  of  the  major  species  (C02, 02,  CO,  H20,  H2,  N2)  based  on 
the  raw  data  (Figures  7  and  8).  Then  we  subtract  a  component  of  the  Raman  signal  equal  to  the  pro¬ 
duct  of  the  correction  factor  for  species  “i”  and  the  value  of  FI: 
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Figure  6.  Scheisati'.’  showing  the  location  of  the  two  additional  photomulti¬ 
plier  tubes  (FI  and  F2)  used  to  monitor  the  laser-induced 
fluorescence  relative  to  other  Raman  channels  in  the  wavelength 
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Figure  9.  Comparison  of  reduced  Raman  data  with  laminar  counterflow 
diffusion  flame  calculation  with  correction  factors  for  fluores¬ 
cence  set  to  zero:  Temperature. 
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Figure  10.  Comparison  of  reduced  Raman  data  with  laminar  counterflow 
diffusion  flame  calculation  with  correction  factors  for  fluores¬ 
cence  set  to  zero:  N2, 02,  and  CO 
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Figure  11.  Comparison  of  reduced  Raman  data  with  laminar  counterflow 
diffusion  flame  calculation  with  correction  factors  for  flores¬ 
cence  set  to  zero:  C02,  H2,  and  CH4. 


(RSc)i  =  (RSuc)i-(CF)iX(FS)Fl 


(2) 


where  RSC  refers  to  the  corrected  Raman  signal,  RSUC  refers  to  the  uncorrected  Raman  signal,  CF  is 
the  correction  factor,  and  FS  refers  to  the  Fluorescence  signal  measured  in  photomultiplier  tube  FI. 
Using  the  estimated  correction  factors  and  Eq.  (2),  the  calculated  temperature  and  mole  fraction 
profiles  are  then  compared  with  the  predicted  values.  The  process  is  iterated  to  convergence.  Fig¬ 
ures  12  to  14  show  the  end  result  of  one  such  process  where  the  inclusion  of  correction  factors  have 
resulted  in  data  agreeing  fairly  well  with  the  predicted  laminar  flamelet  calculations.  An  implicit 
assumption  in  this  calibration  procedure  is  that  the  flames  are  laminar  and  have  very  limited  amount 
of  stretch  such  that  they  do  not  have  any  nonequilibrium  effects.  Predicted  values  for  the  CH4- 
containing  flames  were  obtained  for  different  assumed  stretch  rates  [Chen  et  al.  (1989)].  The  profiles 
shown  in  Figures  7  to  14  were  obtained  for  an  assumed  stretch  of  5/s  and  did  not  depend  strongly  on 
this  assumed  value. 

Once  the  correction  factors  were  obtained  they  were  kept  constant  for  reduction  of  all  the  tur¬ 
bulent  flame  data.  More  sophisticated  methods  using  linear  matrix  manipulation  and  artificial  intelli¬ 
gence  codes  are  being  pursued  elsewhere  pibble  et  al.  (1990)].  Experience  has  shown  that  this  itera¬ 
tion  process  can  be  accelerated  by  “educated  guesses.”  The  calibration  procedures  were  repeated 
before  and  after  each  set  of  measurements. 

23  Discussion  of  Data 

(a)  CO/H2/N2  Flame 

The  fuel  c  rnsists  of  a  mixture  of  27.5%  CO/  32.3%  H2/  40.2%  N2,  a  medium  Btu  gas.  The  exper¬ 
imental  results  are  compared  with  calculations  in  Section  2.5;  here  only  qualitative  observations  are 
made.  The  data  showed  significant  nonequilibrium  chemistry  effects.  However,  even  in  the  high 
shear  region  in  the  neck  of  the  flame  no  significant  localized  extinction  was  observed  because  of  the 
high  hydrogen  content  of  the  fuel.  To  study  a  flame  with  significant  extent  of  local  extinction  it  was 
decided  to  switch  to  pure  methane,  whose  slower  kinetics  are  expected  to  increase  the  probability  of 
localized  extinction.  The  results  of  Masri  et.  al.  (1987, 1988, 1989)  and  Dibble  et  al.  (1987)  support 
this  expectation,  but  Bilger  (1989)  points  out  that  the  reason  may  be  chemical  rather  than  aero¬ 
dynamic. 

(b)  Pure  Methane  Flame 

The  CH4  flame  used  in  this  study  (Figure  5)  is  at  a  Reynolds  number  of  12,000,  based  on  the  jet 
diameter  of  3.175  mm  and  jet  exit  velocity  of  62.5  m/s.  Radial  profiles  of  temperature  and  mole  frac¬ 
tions  of  major  species  were  measured  at  axial  locations  ofx/d  =  5, 10, 20, 30, 50  and  70  and  at  the 
centerline.  The  data  presented  here  were  all  corrected  using  the  factors  obtained  by  the  procedure  of 
Section  2.2.2.  To  help  visualize  the  results  better,  a  schematic  of  bluff-body  stabilized  flame  outlining 
the  various  regions  of  interest  is  shown  in  Figure  15. 

Figures  16  to  19  present  scattergrams  of  temperature  vs.  mixture  fraction  in  various  regions  of  the 
flame.  In  the  recirculation  zone  (x/d  =  5),  most  of  the  data  in  the  T-f  space  lies  on  the  rich  side  close 
to  the  profile  obtained  with  CH4  flames  at  assumed  stretch  of  5/s.  As  Figure  17  shov.s,  the  predicted 
temperature  profiles  for  stretch  of  5/s  and  300/s  are  quite  close  to  each  other.  In  the  neck  region 
(x/d  =  10  -  20)  shown  in  Figure  17  the  scattergram  is  clearly  bimodal.  A  significant  number  of  points 
correspond  to  rich  equilibrated  eddies,  while  a  significant  number  of  eddies  have  fuel-air  mixture 
fraction  values  close  to  stoichiometric  but  significantly  lower  temperature.  These  eddies  manifest 
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Figure  12.  Comparison  of  fluorescence-corrected  Raman  data,  correspond¬ 
ing  to  Fig.  9,  with  laminar  counterflow  diffusion  flame  calcula¬ 
tion. 
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Figure  13.  Comparison  of  fluorescence-corrected  Raman  data,  correspond 
ing  to  Fig.  10,  with  laminar  counterflow  diffusion  flame  calcula¬ 
tion. 
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Figure  14.  Comparison  of  fluorescence-corrected  Raman  data,  correspond¬ 
ing  to  Fig.  11,  with  laminar  counterflow  diffusion  flame  calcula¬ 
tion. 
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Figure  15.  Schematic  of  various  regions  of  the  bluff-body  stabilized 
methane  flame  indicating  various  regions  of  interest 
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Figure  16.  Scattergram  of  temperature  vs.  mixture  fraction  at  x/d  =  5. 


Figure  17.  Scattergram  of  temperature  vs.  mixture  fraction  in  the  narrow 
neck  region  of  the  flame  (x/d  =  10  and  20). 
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Figure  18.  Scattergram  of  temperature  vs.  mixture  fraction  at  x/d  =  30. 


FLUORESCENCE  CORRECTION  FACTORS  INCLUDED 


MIXTURE  FRACTION 


Figure  19.  Scattergram  of  temperature  vs.  mixture  fraction  at  x/d  =  70. 
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localized  extinction  caused  by  high  shear  in  the  neck  region  of  the  flame  (where  the  recirculation  zone 
ends)  and  the  slow  kinetics  of  methane.  Whether  these  eddies  represent  fuel-air  mixture  that  became 
extinguished  because  of  high  strain  or  eddies  that  were  never  ignited  is  difficult  to  determine  based  on 
this  data.  Joint  measurements  of  Raman  and  radical  species  as  OH  are  needed  to  clarify  this  issue. 

At  x/d  =  30  (Figure  18),  there  are  still  a  number  of  locally  extinguished  eddies,  but  the  fuel-air  mix¬ 
ture  is  leaner  and  there  are  some  eddies  which  have  re-ignited  (on  the  lean  side)  and  have  attained 
temperatures  close  to  those  of  strained  flamelets.  As  more  eddies  reach  the  low  shear  region  down¬ 
stream  of  the  neck,  they  appear  to  re-ignite  and  equilibrate  to  high  temperatures,  as  shown  by  data  on 
the  lean  side  in  Figures  18  and  19.  There  are  some  extinguished  eddies  at  x/d  =  18  as  well,  but  then- 
number  is  reduced  significantly.  The  physical  picture  that  emerges  from  the  data  of  Figures  16-19  is 
that  of  a  recirculation  zone  (which  keeps  the  flame  anchored)  connected  to  a  jet-like  region  by  a  fairly 
long  but  narrow  neck  region  of  high  shear  with  a  significant  number  of  locally  extinguished  eddies.  In 
the  neck  region  the  flame  is  intermittently  connected  at  these  high  Reynold  numbers. 

The  scattergrams  of  mass  fractions  of  major  species  in  mixture  fraction  space  in  the  neck  region  of 
the  flame  is  shown  in  Figures  20  to  24  for  02,  N2,  CH*,  H20,  and  C02,  respectively.  Laminar- 
opposed-flow  flame  calculations  for  a  of  5/s  and  300/s  are  also  plotted  in  the  figure  for  reference. 

All  data  has  been  included  in  all  the  scattergrams.  The  results  for  02,  CH4  and  N2  distinctly  show  bi- 
modal  behavior  between  frozen  values  and  equilibrated  values,  whereas  the  data  for  C02  is 
widespread  possibly  because  of  contamination.  It  is  to  be  noted  that  C02  and  CO  species  are  most 
difficult  to  measure  because  of  their  relatively  low  signal  values. 

Finally,  Figures  25  to  29  show  the  measured  profile  of  temperature  and  major  species  across  the 
burner  at  x/d  =  10.  The  temperature  has  been  calculated  based  on  two  methods,  viz.,  sum  of  mole 
fractions  and  Rayleigh  scattering.  The  Rayleigh  signal  is  not  contaminated  by  laser-induced  fluores¬ 
cence  and  therefore  provides  an  independent  temperature  measurement.  Figures  25  and  26  show  that 
the  temperature  measurements  obtained  by  the  two  methods  agree  fairly  well  for  both  mean  and  rms 
values.  The  mean  temperature  peaks  in  the  recirculation  zone  and  has  a  minimum  in  the  core  jet. 

The  rms  temperature  values  peak  in  the  shear  layer  between  the  recirculation  zone  and  the  co-flowing 
air.  The  profiles  for  major  species,  Figures  27-29,  show  reasonable  trends.  There  is  some  entrain¬ 
ment  of  air  in  the  core  of  the  main  jet  though  the  oxygen  is  depleted  there.  Profiles  of  CO  and  H2, 
which  are  intermediates,  peak  in  the  recirculation  zone.  The  mean  values  are  very  low.  It  should  be 
recalled  that  x/d  =  10  corresponds  to  the  narrow  neck  region  containing  a  significant  number  of  high¬ 
ly  strained  eddies.  Finally,  the  products  C02,  and  H20  both  peak  in  the  recirculation  zone  and  are 
fairly  flat  there.  Detailed  comparison  of  data  with  model  predictions  in  the  CH4  flame  is  in  progress. 

2.4  Governing  Equations  and  Models 

This  subsection  describes  the  physical  submodels  and  the  computational  algorithm  used  to  predict 
the  bluff-body  stabilized  CO/H2/N2  flame. 

2.4.1  Navier-Stokes  Equations 

The  governing  equations  for  the  flow  of  a  multicomponent  reacting  gas  are 
•  Continuity: 


&  +  -§-  0*,)  =  0 

dt  dxi  v  ' 


(3) 
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Figure  20.  Scattergram  of  mass-fraction  of  02  vs.  mixture  fraction  in  the 
narrow  neck  region  of  the  flame  (x/d  =  20).  Profiles  obtained 
from  laminar  flamelet  calculations  for  a  of  5/s  and  300/s  are 
also  shown. 
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Figure  21.  Scattergram  of  mass-fraction  of  N2  vs.  mixture  fraction  in  the 
narrow  neck  region  of  the  flame  (x/d  =  10  and  20).  Profiles 
obtained  from  laminar  flamelet  calculations  for  a  of  5/s  and 
300/s  are  also  shown. 
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Figure  22.  Scattergram  of  mass-fraction  of  CH4  vs.  mixture  fraction  in  the 
narrow  neck  region  of  the  flame  (x/d  =  10  and  20).  Profiles 
obtained  from  laminar  flame  let  calculations  for  a  of  5/s  and 
300/s  are  also  shown. 


Figure  23.  Scattergram  of  mass-fraction  of  H20  vs.  mixture  fraction  in  the 
narrow  neck  region  of  the  flame  (x/d  =  10  and  20).  Profiles 
obtained  from  laminar  flamelet  calculations  for  a  of  5/s  and 
300/s  are  also  shown. 
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gure  24.  Scattergram  of  mass-fraction  of  C02  vs.  mixture  fraction  in  the 
narrow  neck  region  of  the  flame  (x/d  =  20).  Profiles  obtained 
from  laminar  flamelet  calculations  for  a  of  5/s  and  300/s  are 
also  shown. 
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Figure  25.  Profile  of  mean  temperature  across  the  burner  at  x/ d  =  10. 

Data  calculated  using  sum  of  mole  fraction  method  and  Ray¬ 
leigh  scattering  are  also  shown. 
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•  Momentum: 


8  ,  s  L  d  ,  s  dp  ^  d  r  ,dui  dui  2  *  dUk . 

-(«)  +  —  («U;)  =  -  ^  +  571^  +  *7  - 


(4) 


Species: 


37c 


dx 


dxj  Sc  dxj 


(5) 


Enthalpy: 


3,,.  3,,v  dp  .  d  r  n  dh  pl  (  Pr  n  ^  ^  ■, 

a  <'*)  +  *T W =  a  +  IT1  Tv  *7 +  Tv (  &  ' 1} "ST1 


dt  dxj  L  Pr  3t;  Pr  v  Sc  '  s%  s  dxj 


(6) 


These  equations  are  supplemented  by  the  equation  of  state 

N 

e-plWY,  (X,/Ms)) 


s  =1 


(7a) 


and  the  thermodynamic  data  relating  the  enthalpy  of  the  mixture  and  of  each  species  to  temperature 

N 


*  =  E  YMT ) 

s=l 


(7b) 


in  which  the  species  enthalpy  hs(T)  is  known. 

Turbulent  flows  are  characterized  by  a  widely  disparate  set  of  length  and  time  scales  so  that  direct 
numerical  solution  woui  i  require  a  very  fine  discretization  of  space  and  time.  Over  a  domain  of 
length  L,  resolution  down  to  the  viscous  length  would  require  on  the  order  of  L/l^  or  Re  number 
of  mesh  points  in  each  direction;  a  similar  argument  holds  for  temporal  resolution.  Introduction  of 
chemical  time-scales  further  exacerbates  the  problem.  Direct  simulations  are  therefore  not  currently 
possible  for  Reynolds  and  Damkol  Jer  numbers  of  interest,  and  indeed  may  never  be.  Turbulence  and 
turbulence-chemistry  interaction  models  are  therefore  required. 

2.4.2  Averaged  Equations  and  Closure 

The  most  common  computational  approach  in  multidimensional  turbulent  flow  is  a  conventional 
or  density-weighted  decomposition  of  the  dependent  variables  followed  by  a  closure  model  for  the 
correlations.  Since  the  Reynolds  number  is  high  in  typical  flows  of  interest  here,  conventional  clo¬ 
sures  assume  that  molecular  transport  processes  are  overwhelmed  by  turbulent  transport.  For  this 
reason,  the  molecular  transport  terms  are  dropped  in  Eqs.  (3)  to  (6),  which  may  eventually  be  written 
as  follows: 
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•  Momentum: 


£  «> +  £  <**> 


dx i  dxj 


(8a) 


Species: 


f  «)  +  £  (ftM 


(8b) 


Here  4  (s  =  1,  Af  +  1)  represents  ail  N  species  mass  fractions  and  the  total  chemical-plus-sensible 
enthalpy. 

The  standard  approach  does  not  account  for  molecular  mixing  processes.  Experiments  in  liquids, 
where  the  Schmidt  number  Sc  is  on  the  order  of  103,  have  indicated  that  molecular  processes  cannot 
really  be  ignored  in  mixing  [Koochesfahani  et  al.  (1985)].  Such  processes  are  less  important  in  tur¬ 
bulent  mixing  of  gaseous  species. 

Closure  is  most  often  accomplished  by  the  k-e  turbulence  model  [Launder  and  Spalding  (1974)] 
with  a  wall-function  treatment  of  the  near-wall  region.  In  the  k-e  model,  the  turbulent  eddy  viscosity 
ft  is  obtained  on  dimensional  grounds  from 

ft  =  C^pk2/  e  (9) 


where  the  turbulence  kinetic  energy  ( k )  and  the  dissipation  rate  (e)  are  obtained  from  modeled  tran¬ 
sport  equations 


_ dk  d  r/ft  .  N  dk  ,  d“i  ft  dp  dp  — 

PUj  ax,  ~  dx,  ^  dx,  ^  ,U  1  dx,  -f  dx,  dx. 


(10) 


and 


de  d  r ,  to  v  de  ,  r  e  ,,  df/,-  ft  J5/T  dp^  —  £1  /ii\ 


dXj  p  dx,  ax. 


The  constants  in  these  equations  are  taken  as  [Launder  and  Spalding  (1974)] 

Cp  =  0.09,  Ci  =  1.44,  C2  =  1.92; 


(12) 


with  turbulent  Schmidt  numbers  for  k  and  e 

ok  =  1.0,  <xe  =  1.30,  (13) 

and  for  thermochemical  scalar  quantities 

a,  =  0.7  (14) 
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The  two-equation  k-e  model  might  be  expected  to  work  in  flows  where  the  Reynolds  stresses  can 
indeed  be  related  solely  to  the  local  strain  rate  [Rodi  (1980)].  If  this  is  not  so,  it  is  necessary  to  ac¬ 
count  for  convection,  diffusion,  and  production/destruction  of  these  stresses,  e.g.,  Amano  and  Kodali 
( 1984).  A  major  advantage  of  jet  flows  is  that  the  k-e  class  of  models  is  more  successful  there  than  in 
recirculating  flows.  The  k-e  model  is  known  to  have  deficiencies  in  flows  with  strong  curvature 
[(Launder  et  al.  (1984);  Gibson,  (1985)]  and  in  strongly  swirling  flows  [(Launder  et  al.  (1984)]  but  is 
successful  in  jet-like  flow  at  least  as  far  as  prediction  of  the  mean  profiles  is  concerned. 

Despite  reservations,  the  k-e  model  is  used  here.  The  long-term  plan  is  to  switch  to  the  Monte- 
Carlo  pdf  transport  model,  which  is  closed  in  the  sense  of  turbulence  correlations. 

We  note  that  all  models  -  k-e,  second-order  Reynolds  stress/flux,  and  pdf  transport  -  obtain 
length-scale  information  from  the  dissipation  rate  equation. 

To  compute  the  transport  of  a  scalar  in  turbulent  flow,  it  is  necessary  to  calculate  the  rate  at  which 
the  fluctuations  in  the  scalar  field  decay.  Most  closure  models  assume  that  the  time  scale  of  scalar 
turbulence  is  proportional  to  that  of  the  turbulent  velocity  field,  usually  characterized  as  qz  j 'e,  where 
q 2  is  twice  the  kinematic  kinetic  energy  of  turbulence  and  e  is  the  energy  dissipation  rate.  Thatmeans 
the  ratio  between  the  velocity  and  the  thermal  time  scales,  R  is  a  constant.  Here  R  =  (g2/e)/(^2/e^), 
where  <f>  designates  a  scalar  variable. 

2.4.3  Chemistry  Model 

The  thermochemistry  must  be  represented  by  as  few  variables  as  possible,  so  as  to  minimize  the 
number  of  field  equations.  Unfortunately,  even  for  the  simplest  hydrocarbons,  ~  50  species  and  ~ 

250  reactions  are  involved.  While  these  detailed  chemical  models  can  be  run  for  one-dimensional 
laminar  flames  (either  premixed  or  diffusion),  their  applicability  to  multidimensional  flames  is  limited 
by  computational  size  and  cost.  Application  to  turbulent  flames  places  further  restrictions  on  the 
schemes.  Thus,  there  is  a  need  for  simplified  yet  reliable  kinetic  models. 

The  partial-equilibrium  model  is  used  here  and  is  presented  in  the  hierarchy  which  starts  with  fast 
chemistry. 

Fast  Chemistry: 

The  scalar  set  <j>s  of  species  mass  fractions  and  mixture  enthalpy  can  be  reduced  to  a  single  con¬ 
served  variable,  £,  called  the  “mixture  fraction”  [Bilger,  in  Libby  and  Williams  (1980)],  as  follows. 
Defining  the  elemental  mass  fraction  z„  which  is  related  to  the  species  mass  fraction  Y}  by 

N 

Zi  =  S  Hj  Yj  1  =  1. M  (15) 

/= i 


where  is  the  mass  of  the  i,h  element  per  unit  mass  of  the  jth  species  in  a  system  of  M  elements  and 
N  species,  and  assuming  equal  diffusion  coefficients,  the  transport  equation  forz,  is 


(16) 


These  elemental  transport  equations  are  homogeneous  since  elements  are  neither  created  nor  des¬ 
troyed  in  chemical  reactions.  Thus  M  transport  equations  result.  With  equal  transport  coefficients,  z, 
is  also  the  total  (chemical  plus  sensible)  enthalpy  which  is  conserved  for  the  low  Mach  number  flows 
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considered  here;  at  high  Mach  numbers  the  mechanical  energy  u2/ 2  cannot  be  neglected,  and  a 
separate  equation  is  needed  for  the  enthalpy;  likewise  for  nonadiabatic  flames. 

If  further,  as  in  the  flows  of  interest  here,  there  are  two  chemically  distinct  reactant  streams  (su¬ 
perscripts  f  for  fuel  and  a  for  air),  the  mixture  fraction  £  is  defined  as 

Z;  —  Z? 

£  =  -7 -  for  each  /  (17) 

z\  -  z  i 


and  vaiies  between  0  and  1  throughout  the  flow. 


The  equation  for  £  is  then  obtained  from  Eqs.  (16)  and  (17). 


and  the  boundary  conditions  are 


£  =  0  in  stream  a 


(18) 


(19) 


£  =  1  in  stream/  (20) 

This  is  a  veiy  popular  approach  [Jones  and  Whitelaw  (1984)]  because  all  of  the  scalar  information  is 
contained  in  the  mixture  fraction  and  can  be  recovered  from  Eq.  (17)  as 

Z,=2?  +  (zf-zf)5  (21) 

and  therefore  only  one  transport  equation,  Eq.  (18),  has  to  be  solved. 

The  atomic  composition  and  total  chemical  plus  sensible  enthalpy  is  sufficient  to  determine  the 
equilibrium  composition  of  the  gas,  by  minimization  of  the  Gibbs  free  energy,  for  constant  pressure 
systems. 

Partial-Equilibrium  Model: 

The  partial-equilibrium  model  adds  one  more  variable  (besides  £)  to  track  the  progress  of  a  pool 
of  intermediate  species  as  the  system  moves  towards  equilibrium.  This  approach  has  been  developed 
for  H2-air  [Janicka  and  Koliman  (1979, 1982);  Bilger  (1980)],  and  CO/H2/N2-air  systems  [Correa  et 
al.  (1984)].  Two-variable  schemes  for  H2-air  and  CO/H2/N2-air  systems  may  be  described  together. 
The  kinetic  mechanism  adopted  consists  of  fast  shuffle  reactions,  which  are  considered  to  be  fast 
enough  for  the  participating  species  to  be  in  partial  equilibrium 

(51)  H  +  02  £  OH  +  O 

(52)  O  +  H2  £  OH  +  H 

(53)  H2  +  0H£H20  +  H 

(54)  20H£H20  +  O 
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(S5)  CO  +  OH^CO,  +  H 


and  relatively  slow  three-body  recombination  reactions 


(Rl)  H  +  OH  +  M  H20  +  M 
(R2)  H  +  O  +  M-OH  +  M 
(R3)  H  +  H  +  M->H2  +  M 
(R4)  CO  +  O  +  M  +  ->  C02  +  M 


Since  the  recombination  reactions  which  govern  the  decay  of  the  radical  pool  are  slow,  it  is  possible 
for  the  radicals  to  exist  in  concentrations  larger  than  when  in  equilibrium.  This  phenome'  ">n  is  re¬ 
ferred  to  as  radical  overshoot  or  superequilibrium  [Correa  et  al.  (1984),  Drake  et  al.  (198/;].  H02 
can  be  included  in  the  scheme  without  difficulty. 

The  kinetic  mechanism  adopted  here  is  simple  so  that  the  turbulence-chemistry  interactions  are 
not  computationally  prohibitive.  It  attempts  to  account  for  radical  overshoot  and  post-flame  burnout. 
Janicka  and  Kollman  (1979)  used  the  scheme  without  reactions  (S5)  and  (R4)  in  their  studies  of  tur¬ 
bulent  hydrogen  flames.  Warnatz  (1981)  demonstrated  its  validity  in  laminar  premixed  hydrogen 
flames  for  temperatures  above  1500  K.  In  earlier  studies,  Kaskan  (1959)  and  Fenimore  and  Moore 
(1974)  used  the  same  shuffle  reactions  to  describe  the  combustion  of  carbon  monoxide  in  laminar 
flames.  Bilger  and  Starner  (1983)  also  discussed  the  merits  of  this  kinetic  scheme  and  concluded  that 
it  would  be  applicable  at  the  higher  temperatures.  At  high  pressures  (well  above  10  atm),  the  rates  of 
the  shuffle  and  recombination  reactions  approach  each  other  and  the  scheme  ceases  to  be  useful. 


To  reduce  the  number  of  conservation  equations  to  be  considered,  a  linear  combination  of  con¬ 
centrations  is  used.  For  this  system,  this  leads  to  the  combined  variable  Y]j2 


Y'h,  -  Y«,  +  T 


M"i  v  ^  MH2 

yoh  + 


2  Mqh 


M0 


♦I 


MHl 

Mh 


Yh  + 


M, 


rco 


CO 


(22) 


The  formation  rate  of  this  variable  becomes  independent  of  the  shuffle  reactions: 

m 

Wff2  =  -  2  MHz  ( wri  +  wR2  +  wR3 ) 


(23) 


Eleven  quantities  are  required  to  describe  the  instantaneous  state  of  the  system,  nine  species  con¬ 
centrations  (02,  Hi,  N2,  CO,  CO i,  Hi0,0,H,0H),  the  density  and  the  temperature.  These  can  all 
be  related  to  just  two  vriables  as  shown  next. 


The  variable  Y//2  varies  between  its  frozen  and  equilibrium  values  leading  to  the  definition  of  the 
reaction  progress  variable  ry. 


Y]{1  -  Yfi2 


n  = 


n2 


(24) 


Hi 


Specification  of  £  amounts  to  five  equations  from  the  definition  Eq.  (17),  while  four  more  indepen¬ 
dent  equations  are  obtained  from  partial  equilibrium  of  the  radical  pool;  specification  of  t?  and  the 
equation  of  state  complete  the  set  of  nonlinear  algebraic  equations  for  the  description  of  the  instan¬ 
taneous  thermochemical  state.  Solution  of  these  equations  gives  a  table  of  concentrations, 
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temperature,  and  density  for  each  £  and  rj.  Representative  results  are  shown  in  Figure  30.  The  max¬ 
imum  temperature  of  ~2170  K  occurs  at  the  stoichiometric  mixture  fraction  0.33,  in  equilibrium  r? 
=  1.0.  The  temperature  decreases  as  rj  decreases  below  the  equilibrium  value  of  unity,  while  OH 
shows  the  expected  superequilibrium  peak.  These  results  are  stored  in  digitized  form  for  use  in  flow 
calculations. 


The  transport  equations  for  the  two  scalar  variables  are 


pui 


,  aUA(rf  JL) 

dxi  dx-t  '  ‘  dxi 


(25) 


pUi 


Mjh 

dxi 


dxi 


)  +  *Ha 


(26) 


showing  the  passive  and  reactive  natures  of  £  and  r?,  respectively.  Equations  (25)  and  (26)  are  analo¬ 
gous  to  the  Navier-Stokes  in  that  they  must  be  subjected  to  time-averaging  and  closure  (or  other 
modeling)  before  becoming  useful  in  the  context  of  turbulent  flames. 

Hydrocarbon  Kinetics: 

Chemical  kinetic  models  are  available  for  the  simple  hydrocarbons  such  as  lower  alkanes,  but 
involve  ~50  species  and  -250  elementary  reactions  if  nitrogen  chemistry  is  included.  Such  schemes 
cannot  be  implemented  in  turbulent  flow  models.  There  is  an  incentive  to  simplify  schemes  while 
preserving  essential  features. 

By  adding  and  subtracting  elementary  reactions  in  the  full  scheme  so  as  to  cancel  out  selected  fast 
reactions,  systematically  reduced  schemes  have  been  obtained  for  simple  hydrocarbons  [e.g.,  Paczko 
et  al.  (1986)].  Typically  three  or  four  progress  variables  are  required,  but  account  for  -50  kinetic 
source  terms  which  are  sums  of  the  surviving  elementary  rates.  Assessment  has  consisted  primarily  of 
running  the  reduced  schemes,  in  calculations  of  laminar  premixed  and  counterflow  diffusion  flames, 
against  the  full  schemes.  In  general,  results  such  as  the  flame  speed  of  premixed  flames  and  the  struc¬ 
ture  of  diffusion  flames  compare  well. 

An  example  of  such  a  scheme  is 


CH4  +  2H  +  H20  ->COt  4H2 

(HI) 

CO  +  h2o  ^  co2  +  h2 

(H2) 

2H2  +  02  -f  2 H20 

(H3) 

3 H2  +  02  ^  2 H20  +  2 H 

(H4) 

Such  models  are  computationally  tractable  in  multidimensional  turbulent  flow  calculations  and  will 
be  used  in  the  future  for  methane  flames. 


32 


2.4.4  Turbulence-Chemistry  Interaction  Models 

Just  as  closure  is  needed  for  velocity-velocity  and  velocity-scalar  correlations,  it  is  an  important  is¬ 
sue  in  the  mean  chenrcal  source  terms.  The  most  widely  used  approach  is  to  derive,  in  some  manner 
to  be  discussed,  the  probability  density  function  for  the  flow  in  terms  of  the  required  scalar(s).  If 
those  scalars  are  (i  =  1,  N),  then  the  average  of  any  quantity  <f>s  is  obtained  by  integration 

&  =  /  •  •  *  /&(«/>; )  P  (x;  rk)  drpi  •••  dtpN  (27) 

where  P  (x;  i/>,)  is  the  pdf  at  the  points  and  <&($)  is  the  instantaneous  value  of  <f>s  corresponding  to 
the  particular  values  (tf>j).  Higher  moments  can  also  be  obtained,  e.g., 

/'?=/ •••  I(<t>s-k)2P&'l’i)drPi  ■■■  difr  (28) 

The  mean  reaction  rate,  mean  density,  or  any  other  scalar  statistic  are  so  obtained.  Also  because 

-  4  p(tfc)  - 

P^^)  =  zrr~P(x-Ad  (29) 

either  density-weighted  or  conventional  statistics  can  be  obtained  for  all  the  <f>s. 

The  closure  problem  is  approached  in  ensemble-averaged  Navier-Stokes  models  by  solving  for  the 
lowest  moments  of  the  pdf  over  scalars.  To  be  specific,  closure  for  the  (two-scalar)  partial- 
equilibrium  model  is  reviewed  next. 

Assumed  Shape  PDF  Model: 

The  mean  and  variance  of  £  and  t?  which  are  obtained  at  each  point  constitute  the  local  (lowest) 
moments  of  the  joint  pdf  of  fluctuations  in  £  and  rj,  Pfcrj).  Moment  methods  assume  that  the  func¬ 
tional  relationship  between  the  joint  pdf  and  these  moments  is  known,  so  that  these  moments  lead  to 
the  local  pdf.  The  functional  relationship  is  assumed  to  be 

(so) 


which  implies  that  fluctuations  in  £  and  rj  are  uncorrelated.  The  pdf  over  £  is  taken  as  a  /?-f  mcdon: 


e q-a*-* 

!<r'  (i-e)1"1^ 

o 


(31) 


In  the  other  dimension  (rj),  Janicka  and  Kollman  (1982)  proposed  the  form  of  the  pdf  as  three 
Dirac  delta  functions: 


HW)  *  Ci  S( i))  +  c2S(ri-fi )  +  c3  fl(l-rj) 


(32) 


The  distribution  consists  of  delta  functions  at  »7=0,  r)=ri  and  r}=  1.  The  constants  a,  b,  and  Cj  -c3  are 
obtained  in  terms  of  the  moments  from  normalization  conditions  and  expressions  for  the  moments 
such  as 

l 

i  =  i  mm  (33) 

0 

where  the  functional  form  of  the  pdf  is  used.  Thus,  the  local  pdf  is  obtained  parametrically  through 
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the  coefficients  which  depend  on  the  local  moments.  These  moment  equations  are  obtained  by 
averaging  the  equations  for  mixture  fraction  and  reaction  progress  variable  [Jones  and  Whitelaw 
(1984),  Correa  and  Shyy  (1987)].  All  mean  quantities  in  the  transport  equations,  e.g.,  wHl,  are  then 
obtained  by  convolution  of  the  instantaneous  quantity,  e.g.,  wHl  (£,r?),  with  the  joint  pdf. 

A  survey  of  the  constants  c  i  -c3  during  calculation  of  jet  flames  has  shown  that  c2  is  very  close  to 
unity.  This  agrees  with  Janicka  and  Kollman’s  (1982)  finding  thatg(r?)  could  be  specified  as  three  del¬ 
ta  functions,  one  delta  function  at  r)=f)  or  a  /?- function  without  much  impact  on  the  predicted  mean 
OH  concentration.  The  adequacy  of  the  shape  of  the  pdf  in  r?  was  tested  by  the  agreement  on  both  the 
mean  and  rms  components  of  the  OH  concentration  [Correa  et  al.  (1984)],  and  in  comparisons  with 
Monte  Carlo  methods  for  direct  calculation  of  the  joint  velocity-scalars  pdf  evolution  equation  [Pope 
and  Correa  (1987)]. 

PDF  Transport  Model: 

A  turbulent  flow  is  described  by  the  following  dependent  variables:  velocity  U  and  additional 
scalars  such  as  mixture  fraction  £  and  the  reaction  progress_variable  ry,  introduced  above.  Here  a 
one-point  statistical  description  in  terms  of  the  joint  pdf  of  U  and  the  scalars  (£  and  rj)  is  sought.  If  the 
flow  is  statistically  stationary,  all  one-point  statistics  depend  only  on  the  spatial  coordinates. 

The  joint  velocity-scalars’  pdf  evolution  model  relaxes  many  of  the  assumptions  made  in  theabove 
standard  closure.  The  joint  pdf  P(V,l,r}\x,r)  is  the  probability  of  the  simultaneous  events!/  =  K,£  = 

£  and  ry  =  rj  at  that  location  in  space  and  time.  All  one-poLit  statistics  are  recovered  from  this  pdf  be¬ 
cause  the  composition  is  a  known  function  of  £  and  rj.  For  example,  the  mean  of  a  quantity  <f>  is 

<<Kx>r)>—  /  M^)  P  (Ktv:x,r)-  dVdldi.  (34) 

Any  correlations  such  as  higher  moments  or  fluxes  can  be  computed  similarly.  The  problem  is  there¬ 
fore  to  determine  the  joint  pdf  P  at  each  point  in  the  flowfield. 

Pope  (1985)  has  developed  a  transport  equation  for  the  evolution  of  the  joint  pdf.  This  exact 
equation  for  P  is  derived  directly  from  the  Navier-Stokes  equations  and  the  species  conservation 
equations.  Even  for  variable-density  flows  with  complex  chemistry,  convection  (including  all  fluxes), 
chemical  reaction,  and  the  mean  pressure  gradient  do  not  have  to  be  modeled.  Molecular  transport 
and  the  fluctuating  pressure  gradient  require  modeling  for  closure.  The  modeling  used  by  Pope  and 
Correa  (1987)  follows  that  for  a  self-similar  plane  jet.  A  conditional  model  was  used  to  discriminate 
between  the  turbulent  and  the  nonturbulent  fluid  in  jet  flows.  The  mixture  fraction  £  also  is  the  con¬ 
centration  of  nozzle  fluid;  hence  £  =  0  indicates  nonturbulent  ambient  gas  and  £  >  0  indicates  tur¬ 
bulent  gas.  Details  of  the  stochastic  mixing  models,  which  account  for  molecular  transport,  and  the 
stochastic  reorientation  models,  which  account  for  pressure  fluctuations,  were  given  first  by  Pope 
(1985). 

The  pdf  transport  occurs  in  a  multidimensional  space  composed  of  the  independent  variables  in 
the  problem,  viz.,  two  spatial  variables,  ( x,r ),  three  velocity  variables  V,  and  two  scalar  variables  (£,ry). 
Calculations  of  transport  equations  are  already  computationally  intensive  in  three  dimensions  and  so  a 
Monte  Carlo  numerical  technique  is  used  [Pope  (1985)].  Simplifications  analogous  to  the  boundary- 
layer  approximation  were  made. 

In  the  Monte  Carlo  method,  the  joint  pdf  is  represented  by  a  large  number  N  of  notional  particles 
(iV  ~  35,000).  At  location*,  the  «-th  particle  has  the  radial  position  r(”\x),  velocity  U(n\x)  and  com¬ 
position  £("',  rfn\x).  Initial  conditions  are  specified  at  the  jet  exit  (*  =  0). 
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This  model  has  been  applied  to  the  Re =8500, 40%  CO,  30%  H2, 30%  N2  turbulent  jet  diffusion 
flame,  with  partial-equilibrium  chemistry.  Although  the  results  agreed  broadly  with  those  of  the 
assumed-shape  pdf/k-e  model,  some  differences  were  evident  [Pope  and  Correa  (1987)].  For  exam¬ 
ple,  a  strong  correlation  was  found  to  exist  between  £  and  rj,  but  it  did  not  greatly  affect  the  predicted 
mean  quantities.  The  prediction  of  rms  quantities  was  improved. 

2.4 .5  Flow  Algorithm 

This  account  of  the  flow  algorithm  is  brief,  since  it  has  been  well  published  already.  Further  de¬ 
tails  and  applications  may  be  found  in  the  review  by  Correa  and  Shyy  (1987). 

Axisymmetric  bluff-body  stabilized  flames  are  characterized  by  recirculating  rather  than  jet-like 
flow,  with  attendant  numerical  complications.  The  governing  equations  are  written  in  strong  conser¬ 
vation  law  form  in  general  curvilinear  coordinates,  as  derived  by  Vinokur  (1974).  The  conservation 
equations  can  be  written  in  Cartesian  coordinates  for  a  dependent  variable  <f>  in  the  following  form,  re¬ 
taining  3D  formulation  for  the  present: 

ii  w) +  fr  <***) +  i  f : 1  +  f  f : > +  i (r  §: 1 + R  (w)  (35a) 


where  F  is  the  effective  diffusion  coefficient  and  R  is  the  source  term.  When  new  independent  vari¬ 
ables  £,ri,  and  7  are  introduced,  Eq.  (35a)  changes  according  to  the  general  transformation 
C  =  £(x,y,z),  t)  =  rj(x,y,z),  7  =  7 (x,y,z).  The  result  of  this  coordinate  transformation  is  to  transform 
the  arbitrarily  shaped  physical  domain  into  a  rectangular  parallelopiped. 


Equation  (35a)  can  be  rewritten  in  (£,17,7)  coordinates  as  follows: 

w) + ^  m)  *  ^  (m) 


(35b) 


d_ 

% 


J  (<?!!  <k  +<I  V2  <f>n  +  <713^7)1 


+  dr, 


y  (?2i  <k  +  922  +  923  ^7) 


+  dl 


J  (?31  <k  +  932  4>n  +  933  <f>t) 


where  the  contravariant  velocity  components  and  the  combinations  of  metric  terms,  q l},  are 
U  =  u(y^  -  y^zn)  +  v(x^zn  -  xfa)  +  w(xy^  -  x^y„) 


V  =  u(y^z(  -ypj  +  K*«z7  -x*i)  +  w(x^(  ~  x$i ) 


(36) 


(37) 


W  =  M(y€z„  -  y „z€)  +  v(x^z€  -  XfZn)  +  -  *^)  (38) 

9 11  =  OV7  "  yfn?  +  (Wn  ~  xt?t)2  +  (-V 7  "  Wn?  (39) 
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? 22  =  -  yfa?  +  (xfa  -  x^z€)2 3 4+  -  X^7)2  (40) 

4  33  =  O^r,  -  y#z)2  +  (V«  “  xtzl)2+  (x&r,  ~  *r^f)2  (41) 

912  -  (y Tf2^  “  y (,z-i)  (yifi  ~  y-rzn)+  ^izt  ~xrz()  ix~rzn  ~  xnzt)  "*"  (xt ~  x&i)  ~  xT?n)  (42) 

4  13  ~  (y«Z-7  -  (y^f  -  y^ >j)  +  -  Xrjzi)  ixrf(  -  •z{zij)+  (*^>1  -  xiy()  (^,  -  x'y’n)  (43) 

<723  =  (y«z,  -  y^)  6v*  -  >-,.2,)+  (x,?(  -  X(zn)  (x^  -  x&)  +  (X(yn  -  xy()  (xj(  -  X(yJ  (44) 

with  q2!  =  q12,  q3i  =  qi3.  q32  =  q23,  and  finally 

J  =x^y  fa  +  xyfa  +  -  x$fn  -  x^fz  -  xtffa  (45) 


S  (f,  27,7)  is  the  source  term  of  the  governing  equation  in  (£,  >7,7)  coordinates.  The  momentum  equa¬ 
tions  are  written  in  terms  of  the  Cartesian  velocity  components,  while  in  the  continuity  equation  the 
contravariant  velocity  components  are  used.  This  formulation  has  the  advantage  that  the  equations 
retain  simple  forms  without  a  large  number  of  additional  terms. 

The  main  features  of  the  computational  method  are 

1.  Meshing:  The  domain  of  interest  in  the  physical  domain  is  covered  with  a  body-fitted  computa¬ 
tional  mesh  and  is  numerically  transformed  to  a  square  in  2D  (cube  in  3D).  The  transformed 
equations  (35b)  are  solved  in  the  transformed  domain.  Of  course,  most  laboratory  experiments, 
including  the  present  one,  can  be  described  in  terms  of  a  rectilinear  (Cartesian  or  cylindrical- 
polar)  mesh,  perhaps  nonuniform  in  cell  size. 

2.  Staggered  Grid:  Variables  are  not  defined  at  the  same  mesh  points.  The  “staggered”  grid  con¬ 
cept  [Harlow  and  Welch  (1965)]  locates  scalar  variables  such  as  pressure,  density,  turbulence 
kinetic  energy,  and  so  on,  at  common  mesh  points.  Offset  half  a  mesh-spacing  to  the  left  and 
below  are  the  (axial)  u  and  (lateral)  v  velocity  points,  respectively  (in  three  dimensions,  w  is 
similarly  staggered).  This  configuration  stabilizes  the  pressure-velocity  coupling. 

3.  Control  Volume  Formulation:  The  partial  differential  equations  are  integrated  over  control 
volumes  centered  on  each  grid  point.  The  advantages  of  the  staggered  grid  are  then  realized, 
since  the  velocity  and  scalar  values  are  located  exactly  where  they  are  needed  -  at  the  surfaces 
of  the  control  volumes.  This  formulation  also  ensures  global  conservation. 

4.  Discretization  Schemes:  Discretization  of  the  convection  terms  can  be  one  of  the  principal 
sources  of  error  in  numerical  solutions  of  high  Reynolds  number  flow.  Second-order  central 
differences  for  the  convection  term  are  consistent  with  discretization  of  diffusion  terms,  but 
make  the  coefficient  matrix  lose  diagonal  dominance  when  the  cell  Reynolds/Peclet  number 
exceeds  two.  Solution  profiles  then  exhibit  characteristic  “wiggles”  or  oscillations  which  can 
become  unbounded.  A  standard  method  to  stabilize  high  Reynolds/Peclet  number  calculations 
is  to  evaluate  the  fluxes  using  the  “hybrid”  scheme,  consisting  of  first-order  upwind  differencing 
at  cells  where  the  Reynolds/Peclet  number  exceeds  two  and  second-order  central  differencing 
everywhere  else  [Spalding  (1972)].  The  use  of  the  hybrid  scheme  therefore  causes  the  loss  of 
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one  order  of  accuracy  in  regions  of  high  (greater  than  two)  cell  Reynolds  or  Peclet  numbers, 
leading  to  overly  diffusive  solutions  [Leschziner  and  Rodi  (1980);  Raithby  (1976);  Shyy  and 
Correa  (1985)].  Numerical  stability  is,  however,  guaranteed. 

The  first-order  upwinding,  second-order  upwinding,  and  QUICK  [Leonard  (1981)]  schemes 
are  coded  into  our  algorithm,  and  any  one  of  them  can  be  used  in  any  given  transport  equation 
by  selecting  appropriate  input  options. 

Evaluations  of  various  finite  difference  operators  for  the  convection  term  in  the  context  of 
two-  and  three-dimensional  recirculating  flow  calculations  can  be  found  in,  for  example,  Tuann 
and  Olson  (1978),  Leschziner  and  Rodi  (1980),  Shyy  and  Correa  (1985),  and  Mei  and  Plotkin 
(1986).  Discretization  errors  have  also  been  quantitatively  compared  to  other  terms  in  each 
balance  equation  [Correa  (1984),  McGuirk  et  al.  (1981)].  A  second-order  centrally-differenced 
a  posteriori  calculation  of  convection  terms  is  subtracted  from  the  original  operator.  The  out- 
of-balance  term  is  the  “numerical  diffusion”  and  is  compared  with  each  of  the  other  terms  (con¬ 
vection,  diffusion,  pressure  gradient,  source)  at  each  node.  Correa  (1984)  has  already  provided 
such  evaluations  for  the  bluff-body  stabilized  diffusion  flame  of  Roquemore  et  al.  (1982)  and 
has  indicated  where  to  concentrate  the  grid  in  such  a  calculation. 

In  the  calculations  reported  below,  contours  of  the  transverse  cell  Reynolds  number  are 
presented  to  help  quantify  the  errors  due  to  discretization. 

5.  SIMPLE:  In  the  set  of  fluid  flow  equations,  pressure  is  the  only  variable  that  does  not  have  a 
governing  equation.  Furthermore,  at  low  Mach  numbers  the  pressure  is  essentially  constant 
and  so  the  continuity  equation  cannot  be  solved  directly.  The  so-called  “Semi-Implicit  Method 
for  Pressure-Linked  Equations”  (SIMPLE)  algorithm  [Patankar  (1980)]  is  intended  to  over¬ 
come  this  problem.  The  essential  feature  is  the  replacement  of  the  continuity  equation  (which 
does  not  contain  the  pressure)  with  a  pressure  correction  equation,  and  subsequent  sequential 
manipulations  of  the  velocity  field. 

The  algorithm  proceeds  sequentially  through  the  momentum,  pressure-correction,  and  other 
scalar  equations  with  a  maximum  allowance  being  prescribed  for  updating  each  dependent  vari¬ 
able.  This  procedure  is  defined  as  the  outer  iteration  cycle.  Within  each  equation,  and  to 
within  a  prescribed  tolerance,  the  discretized  system  of  linearized  equations  is  also  solved  itera¬ 
tively.  This  procedure  is  called  the  inner  iteration  cycle.  The  degree  of  convergence  of  the  solu¬ 
tion  of  the  individual  equations  in  the  inner  cycle  can  be  very  influential  on  the  overall  rate  of 
convergence  of  the  outer  cycle. 

6.  TDMA:  The  algebraic  equations  that  are  the  result  of  discretization  are  cast  into  a  tridiagonal 
form  and  solved  by  line  relaxation  using  the  tridiagonal  matrix  algorithm  (TDMA). 

2.5  Comparison  of  Model  Predictions  with  Data 
2.5.1  Grid 

The  nonuniform  but  rectilinear  75  (axial)  x  60  (radial)  grid  used  for  the  axisymmetric  2D  calcula¬ 
tion  is  shown  in  Figure  3i.  In  each  of  the  coordinate  directions,  the  mesh  spacing  grows  according  to 
a  geometric  progression  with  a  modest  expansion  i  atio  (2%  in  x-direction,  3%  in  y-direction).  The 
domain  measures  55  jet  diameters  in  axial  extent,  well  beyond  the  high  shear  region  of  interest.  The 
transverse  grid  is  chosen  to  coincide  with  the  fixed  geometric  points  defined  by  the  fuel  tube  radius, 
the  bluff-body  radius,  and  the  tunnel  radius.  Since  the  tunnel  is  actually  of  square  cross  section,  the 
outer  radius  was  defined  to  give  the  same  cross-sectional  area  as  the  square.  The  tunnel  wall  has  no 
direct  effect  on  the  flame  since  the  flame  is  confined  to  the  vicinity  of  the  centerline;  there  is  a  weak 
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Figure  31.  Nonuniform  rectilinear  75  (axial)  x  60  (radial)  grid. 
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effect  through  the  axial  pressure  gradient  caused  by  the  heat  release,  but  this  effect  is  accounted  for  by 
preserving  cross-sectional  area  of  the  duct. 

The  grid  expands  in  the  radial  and  the  downstream  axial  directions,  but  is  dense  enough  in  the 
critical  shear  regions  as  shown  in  subsection  2.5.3. 

252  Boundary  Conditions 

The  burner  conditions  are:  pressure  =  1  atm,  inlet  temperature  =  300  K,  bulk-averaged  jet  veloci¬ 
ty  =  80  m/s,  and  bulk-averaged  annulus  velocity  =  6.5  m/s.  Boundary  conditions  in  the  calculation 
are  taken  as  plug  flow  in  the  jet  and  annulus,  with  inlet  turbulence  kinetic  energy  ( k )  taken  as  0.5%  of 
the  mean  local  velocity  squared,  and  inlet  dissipation  rate  (e)  obtained  from  the  assumed  length  scale 
of  turbulence  (3%  of  the  jet  or  annulus  dimension).  It  is  well  known  that  such  calculations  can  be 
sensitive  to  inlet  boundary  conditions  on  k  and  e,  but  variations  in  a  reasonable  range  made  little 
difference.  For  example,  increasing  the  jet  turbulence  length  scale  by  a  factor  of  5  changed  the  mean 
mixture  fraction  field  values  by  less  than  1%.  The  reason  is  that  turbulence  generation  in  this  flow  ra¬ 
pidly  overwhelms  the  inlet  levels. 

2.53  Reduction  of  Discretization  Error 

As  discussed  in  the  section  on  “Discretization”  above,  numerical  diffusion  errors  are  expected  in 
regions  where  the  cell  Reynolds  numbers  (Peclet  numbers  in  scalar  transport  equations)  are  above 
two.  Figure  32  shows  the  transverse  cell  Reynolds  numbers  in  the  flame  calculation.  The  local 
transverse  cell  Reynolds  number  at  eacn  point  (x,y)  in  the  calculation  domain  is  defined  as 


Rc„  (rv\  -  P(x<y)vMty 

R  Ayi  ,y)  ~  tt(x,y)  v46) 

The  transverse  direction  is  significant  since  it  is  the  direction  in  which  the  important  turbulent 
exchange  (of  mass,  momentum,  energy)  occurs.  The  density,  velocity,  and  turbulent  viscosity  fields 
are  taken  from  the  converged  solution.  It  is  clear  from  Figure  32  that  the  cell  Reynolds  numbers  are 
below  two  everywhere  except  in  the  shear  layer  coming  off  the  bluff  body:  this  is  not  a  cause  for  con¬ 
cern  since  there  is  no  fuel  or  flame  there  (as  shown  below).  The  important  shear  layers  are  those 
along  the  flammable  contours  of  mixture  fraction,  and  these  are  adequately  gridded. 

Eliminating  numerical  diffusion  from  the  convection  operator  does  not  mean  that  there  are  no 
numerical  errors.  A  mesh-size  dependent  error  still  exists  because  of  truncation,  i.e.,  neglecting  terms 
of  higher-order  (than  second)  in  the  Taylor  series  expansion  of  the  convection  and  diffusion  operators 
[Correa  and  Shyy  (1987)].  It  does  mean,  however,  that  the  errors  will  monotonically  decline  in  a 
sequence  of  nested  grids,  and  also  that  turbulence  models  can  be  assessed  without  artificial  diffusion. 

2.5.4  Contour  Plots 

Contour  plots  of  the  calculated  fields  are  presented  first  to  give  an  overall  view  of  the  flame  struc¬ 
ture.  Profiles  are  compared  with  data  in  the  following  two  subsections. 

Figures  33  and  34  show  contours  of  the  Favre-averaged  mean  mixture  fraction  and  its  variance, 
respectively,  which  are  the  lowest  moments  of  the  pdf  over  mixture  fraction.  For  reference,  note  that 
the  stoichiometric  value  is  £ = 0.33.  Both  sets  of  contours  appear  jet-like  but  with  greater  mixing  at  the 
base  ( 0  <  x/d  <  20).  The  maximum  strain  rates  (scalar  fluctuations)  occur  in  the  jet  shear  layer  and 
are  on  the  order  of  30%  of  the  local  mean.  Contours  of  the  turbulent  eddy  viscosity  (&),  obtained 
from  the  k-e  model  as  described  above,  show  that  there  is  another  shear  layer  associated  with  the  bluff 
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Figure  32.  Transverse  cell  Reynolds  numbers  in  the  calculation. 
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Figure  33.  Contours  of  the  Favre-averaged  mean  mixture  fraction. 
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Figure  34.  Contours  of  the  Favre-averaged  mixture  fraction  variance. 


body.  The  latter  layer  affects  only  the  velocity  (and  so  k,  e,  and  &)  since  the  mixture  fraction  is  not 
transported  that  far  out  radially. 

The  peak  mean  (unweighted)  temperatures  (Figure  35)  are  about  500  K  below  the  stoichiometric 
adiabatic  equilibrium  value  of  2170  K  This  large  decrement  is  caused  by  the  intense  turbulence.  The 
reaction  progress  variable  r?  (Figure  36)  likewise  shows  a  minimum  in  the  shear  layer. 

Mean  species  fields  are  shown  in  Figures  37.  The  minor  species  peak  along  the  stoichiometric 
contour,  as  expected. 

2.5.5  Axial  Profiles 

Axial  profiles  along  the  centerline  are  presented  for  the  mean  density,  axial  velocity,  turbulence 
kinetic  energy,  and  the  mean  a  J  variance  of  scalar  fluctuations  (Figures  38  to  42).  Raman  data  are 
shown  for  the  scalar  fields.  The  predicted  jet  decays  faster  than  measured.  This  error  can  be  attri¬ 
buted  primarily  to  use  of  the  standard  k-e  turbulence  model,  since  numerical  diffusion  in  the 
transverse  direction  is  largely  eliminated  as  shown  above.  Deficiencies  of  the  k-e  model  are  well  docu¬ 
mented  in  the  literature  and  include  the  assumptions  of  isotropic  turbulent  eddy  viscosity  and  gradient 
diffusion,  the  need  for  semi-empirical  modification  in  axisymmetric  jet  or  jet-dominated  flow,  and  the 
intrinsic  failure  to  account  for  transport  of  turbulent  stresses.  Overdiffusive  predictions  by  the  k-e 
submodel  seem  to  be  responsible  for  the  predicted  jet  decaying  too  rapidly. 

Figure  43  shows  that  the  normalized  turbulence  intensity  v  (2fc/3)  /  U,  where  U  is  the  local  mean 
axial  velocity,  asymptotes  to  ~0.2.  The  corresponding  normalized  scalar  fluctuation  intensity  profile 
(rms/mean  scalar)  is  compared  with  Raman  data  in  Figure  44.  The  scalar  fluctuations  develop  too  ra¬ 
pidly  in  the  calculation,  being  about  twice  the  measured  value  at  x/d=  15,  but  the  asymptotic  value  of 
0.3  agrees  with  the  data.  It  is  interesting  to  note  that  an  asymptotic  value  of  0.3  has  been  calculated 
for  turbulent  jet  diffusion  flames,  also  in  agreement  with  Raman  data  on  scalar  fluctuations  (Lapp  et 
al.  (1983)]. 

Note  that  second-order  closure  of  the  turbulence  correlations  could  have  been  employed  to  relax 
some  of  the  /c-e-scalar  variance  model  assumptions.  The  plan,  however,  is  to  proceed  directly  to  the 
Monte-Carlo  joint  pdf  transport  method,  which  is  closed  to  all  orders  and  will  provide  the  same  level 
of  modeling  sophistication  in  recirculating  flow  as  already  demonstrated  for  jet  flames  [Pope  and 
Correa  (1987)]. 

2.5.6  Radial  Profiles 

Figure  45  compares  various  Raman  data  with  the  calculations  aix/d=  10,  which  is  in  the  region  of 
maximum  shear.  It  is  again  clear  that  peak  temperatures  are  well  below  the  stoichiometric  adiabatic 
equilibrium  value.  Species  predictions  are  generally  in  reasonable  agreement  with  the  data,  except  in 
the  cool  fuel-rich  core  where  partial  equilibrium  in  the  radical  pool  breaks  down.  There  the  model 
predicts  too  much  H2O  at  the  expense  of  H2.  Breakdown  of  the  partial  equilibrium  thermochemical 
model  has  previously  been  noted  in  our  jet  flame  calculations,  with  both  assumed  shape  as  well  as 
Monte  Carlo  pdf  methods,  and  is  to  be  expected  since  the  assumption  of  partial  equilibrium  in  the 
radical  pool  is  invalid  below  -1100  K. . 

Overprediction  of  the  jet  decay  rate  can  be  interpreted  retrospectively  as  a  shift  in  virtual  origin. 
Figure  46  compares  the  calculation  &tx/d=  15.4  with  the  data  at.r/d=20.  This  axial  shift  results  in  ex¬ 
cellent  agreement  on  the  mean  mixture  fraction  and  the  temperature  fields. 

Figures  47  and  48  compare  the  predicted  Favre-averaged  rms  mixture  fraction  fluctuations,  nor¬ 
malized  by  the  local  mean,  with  the  Raman  data  at x/d=  10  andx/rf=40,  respectively.  This  normaliza- 
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Figure  36.  Contours  of  the  mean  reaction  progress  variable. 
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Figure  37.  Contours  of  the  Favre-averaged  mean  species  mass  fractions, 
(a)  H2. 
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Figure  37.  Contours  of  the  Favre-averaged  mean  species  mass  fractions 

(c)  N2. 

(d)  CO. 
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Figure  37.  Contours  of  the  Favre-averaged  mean  species  mass  fractions. 

(e)  O. 

(f)  OH. 


47 


Figure  37.  Contours  of  the  Favre-averaged  mean  species  mass  fractions. 

(g)  H. 

(h)  H20. 
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Figure  37.  Contours  of  the  Favre-averaged  mean  species  mass  fractions. 

(i)  co2. 


Figure  38.  Centerline  axial  profile  of  the  mean  density. 
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Figure  45.  Comparison  of  predicted  radial  profiles  with  Raman  data  at 
x/d  =  10. 

(a)  Mean  temperature. 

(b)  Mean  mixture  fraction. 
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Figure  46.  Comparison  of  predicted  radial  profiles  at  x/d  =  15.4  with 
Raman  data  at  x/d  =  20. 

(a)  Mean  mixture  fraction. 

(b)  Mean  temperature. 


tion  takes  care  of  the  problem  with  the  virtual  origin  but  towards  the  edges  of  the  mean  mixture  frac¬ 
tion  profile,  where  the  mean  tends  to  zero,  the  quantity  becomes  indeterminate.  The  agreement  at 
x/d=  10  is  fairly  good,  although  the  off-axis  peak  is  overpredicted  by  about  50%  consistent  with  the 
overprediction  of  the  near-exit  turbulence  levels.  The  agreement  at  x/d=40  is  very  good. 

2:6  Monte-Carlo  Pdf  Transport  Method  in  Recirculating  Flow 

The  assumption  of  shape  -  particularly  the  assumption  of  statistical  independence  of  the  scalar 
fluctuations  -  limits  the  number  of  scalars  in  the  assumed  shape  pdf  method.  Direct  pdf  evolution 
methods  are  intrinsically  suited  to  this  problem  of  closing  the  chemistry.  They  can  also  relax  many  of 
the  turbulence  model  assumptions,  such  as  gradient  diffusion,  made  above.  Substantial  difficulties  ar¬ 
ise,  however,  in  recirculating  flow  particularly  in  complex  geometries  because  the  wall  shape  sets  the 
mean  pressure  field. 

A  derivative  of  the  full  pdf  evolution  method  is  to  use  conventional  k-e  or  Reynolds  stress  closure 
for  the  momentum  equations,  and  the  pdf  evolution  method  for  the  scalar(s)’  pdf  only.  Gradient 
diffusion  is  therefore  invoked  in  the  scalar  flux,  but  there  is  no  need  for  assumption  of  pdf  shape.  This 
approach  is  useful  in  internal  recirculating  flows  where  the  mean  pressure  field  is  very  important  in 
setting  the  mean  velocity  field,  as  in  bluff-body  stabilized  burners.  This  mean  pressure  field  is  not 
provided  by  the  pdf  evolution  method.  To  relax  even  the  gradient  diffusion  assumption,  the  full  pdf 
method  could  be  used  with  (i)  the  mean  pressure  field  given  by  the  mean  momentum  equation  and 
(ii)  the  scalar  fluxes  given  by  the  velocity-scalar  coupling  in  the  pdf,  beyond  which  the  velocity  infor¬ 
mation  in  the  pdf  would  be  disregarded.  The  procedure  is  iterated  to  convergence. 

Preliminary  work  on  the  hybrid  approach  has  begun,  with  a  goal  of  reproducing  the  CO/H2/N2 
flame  calculation  discussed  above,  and  then  going  on  to  more  complex  chemistry  (methane).  A  calcu¬ 
lation  procedure  which  iterates  between  the  hydrodynamics  and  the  pdf  evolution  equations  has  been 
worked  out.  The  hydrodynamical  part  contains  the  mean  flow  and  turbulence  model  equations,  and 
the  geometry  of  the  problem.  The  mean  velocity  and  the  turbulence  kinetic  energy  fields  are  passed 
to  the  pdf  evolution  equation  solver.  The  pdf  is  shifted  in  phase  space  to  agree  with  these  lower  mo¬ 
ments. 

The  test  case  is  an  axisymmetric  constant  density  1:2  sudden-expansion  flow.  The  inlet  boundary 
conditions  are  a  plug  100  m/s  axial  velocity  and  a  step-function  profile  for  the  conserved  scalar.  This 
(purely  illustrative)  calculation  used  a  21  x  21  grid  for  the  mean  flow.  Figure  49  shows  contours  of  the 
stream  function,  and  the  mean  and  variance  of  the  conserved  scalar,  all  of  which  behave  as  expected. 

The  next  step  is  to  add  scalar  dimensions  to  the  pdf  and  redo  the  two-scalar  flame  calculation  of 
Section  2.5.  The  flexibility  to  handle  multiscalar  kinetics,  without  realizabililty  constraints  imposed  by 
assumed  shape  pdfs,  will  be  exploited  to  model  the  methane  flame  as  well  (three  or  four  scalars,  as 
necessary). 

2.7  Universality  of  Turbulence-Chemistry  Models 

Should  turbulent  flames  be  viewed  as  ensembles  of  strained  laminar  “flamelets”  [Peters  (1987)]  or 
as  broader  “distributed  zones”  [Bilger  (1989)]  of  reacting  species?  This  section  shows  that  the 
answer  is  not  unique  even  in  a  given  flame.  Turbulence-chemistry  interactions  span  several  orders  of 
magnitude  in  Damkohler  number,  depending  on  the  specific  reaction  and  the  specific  turbulence  pro¬ 
cess  in  question.  Extreme  cases  are 

(i)  Relatively  weak  interactions  which  are  characterized  by  moderate  Reynolds/Peclet  numbers, 
yielding  order-of-unity  Damkohler  numbers  for  the  radical  recombination  reactions  and  large  Dam- 
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Figure  49.  Application  of  hybrid  Monte-Carlo-pdf-transport/finite- 

diJEference  code  to  axisymmetric  constant-density  V2  sudden- 
expansion  flow  with  step-function  inlet  profile  for  the  conserved 
scalar. 

(a)  Stream  function,  contour  spacing  0.03. 

(b)  Mean  of  conserved  scalar,  contour  spacing  0.1. 

(c)  Variance  of  conserved  scalar,  contour  spacing  0.01. 
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kohler  numbers  for  the  chain-branching  reactions.  This  regime  exhibits  superequilibrium  radical  lev¬ 
els,  which  control  burnout  of  major  species  and  CO  as  well  as  oxidation  of  N2. 

(ii)  Strong  interactions,  where  even  ignition  reactions  may  have  order-of-unity  Damkohler 
numbers  leading  to  a  question  of  localized  turbulence-induced  extinction  and  eventually  to  blowoff. 

To  illustrate  the  range  of  interactions,  chemical  and  fluid-mechanical  scales  have  been  compared 
for  two  Tpes  of  air-breathing  combustors.  Diffusion  flames  have  no  intrinsic  length  scale,  so  time 
scales  u  e  used  to  compare  mixing  and  chemistry.  The  physical  picture  to  bear  in  mind  is  whether,  at 
an  interface  between  dissimilar  reactants,  turbulent  straining  occurs  rapidly  or  slowly  relative  to  the 
rates  of  important  reactions. 

Selected  characteristic  kinetic  times  were  obtained  from  stirred  reactor  calculations,  under  condi¬ 
tions  of  (i)  0.1  atm  H2-air  combustion,  1  ms  residence  time,  H2  at  1100  K,  air  at  800  K,  equivalence 
ratio  <j>  =  0.8,  and  (ii)  10  atm  CH^-air  combustion,  10  ms  residence  time,  CH4  at  300  K,  air  at  610  K, 
and  <p  =  0.7,  a  stoichiometry  representative  of  low-NO*  “head-end”  conditions.  The  two  cases  ap¬ 
proximate  a  hypothetical  supersonic  combustion  ramjet  (“scramjet”)  in  the  mid-Mach  number  range 
(10-15)  and  a  methane-fueled  gas  turbine  combustor,  respectively,  closely  enough  for  the  present  pur¬ 
pose  of  extracting  kinetic  time  scales. 

The  kinetic  scheme,  adopted  from  Glarborg  et  al.  (1986),  contains  25  species  and  100  reactions  for 
methane  oxidation  excluding  C2-chemistry.  The  “thermal”  and  “prompt”  (alkyl  species-plus- 
nitrogen)  NO*  mechanisms  are  also  included  because  pollutant  chemistry  introduces  additional  time 
scales.  Concentrations,  density,  and  temperature  from  each  of  the  stirred  reactor  solutions  were  used 
to  estimate  characteristic  times  of  the  ten  reactions  selected  in  Table  I.  Characteristic  chemical  time 
scales  are  estimated  as  rc  =  l./(ky[X])  where  kf  is  the  reaction  rate  and  [X]  is  the  concentration  of  the 
more  abundant  reactant,  with  an  additional  concentration  factor  for  three-body  reactions.  Time 
scales  for  selected  reactions  are  shown  in  Table  I. 

Turbulence  scales  are  estimated  as  follows.  The  integral  time  scale  is  given  by  77  =  A/u'  where  u ' 
is  the  rms  velocity  fluctuation  and  A/  is  the  integral  length-scale.  From  the  kinematic  viscosity  v  and 
the  dissipation  rate  e=u'3/A/,  the  Kolmogorov  time  scale  tk  follows  as  tk  =  (iz/e)1/2.  The  relevance 
of  the  Kolmogorov  time  scale  is  worth  attention.  tk  characterizes  the  intense  turbulent  straining  asso¬ 
ciated  with  the  smallest  scales.  Since  the  pdf  for  scalar  dissipation  is  log-normally  distributed,  this 
(tk)  process  is  highly  intermittent  and  is  experienced  by  only  a  small  part  of  the  flow  at  any  given  in¬ 
stant. 

These  estimates  are  applied  :o  the  two  cases  described  above,  with  the  following  assumptions:  (i) 
in  the  supersonic  combustor,  the  mean  axial  velocity  <u>  =  3000  m/s,  u'  =  300  m/s,  A/  =  0.002  m .; 
the  latter  follows  ffcrn  PLIF  (planar  laser -induced  fluorescence)  flow  visualization  of  the  mixing  of 
jets  injected  behind  a  realistically  sized  step  into  a  supersonic  crossflow  [Correa  et  al.  (3989)]  and  (ii) 
the  mean  velocity  <u>  =  100  rn/s,  u'  =  10  m/s,  A/  =  0.005  m,  being  typical  conditions  in  a  turbine 
combustor.  Turbulence  time  scales  are  shown  in  Table  II. 

Numerical  values  in  Tables  I  and  II  may  be  used  to  compute  the  turbulent  Damkohler  numbers 
jy/rc  and  tk/tc,  the  latter  also  being  an  inverse  Kaiiovitz  number,  and  so  to  compare  the  turbulence 
and  chemical  scales: 

(i)  Scramjet:  Because  of  the  long  chemical  time  scales  (caused  by  the  low  static  pressures),  the 
reactions  are  slow  relative  to  both  the  integral  and  Kolmogorov  scales.  The  exceptions  are  the  two- 
body  oxyhydrogen  reactions  which  are  comparable  to  the  integral  scale.  The  recombination  reaction 
time  scales  are  comparable  to  the  residence  time,  which  is  why  there  is  a  combustion  efficiency  issue 
in  such  devices.  The  ordering  of  time  scales  implies  that  pastial  equilibrium  in  the  radical  pool  can  be 
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a  useful  simplifying  assumption.  In  any  event,  it  is  unlikely  that  the  flamelet  model  would  be  applica¬ 
ble.  Flame  broadening  by  finite-rate  chemistry  would  destroy  the  laminar  flamelet  topology. 

(ii)  Turbine  combustor:  Here  the  oxyhydrogen  shuffle  reactions  are  fast  compared  with  the  tur¬ 
bulence,  except  for  O  +  H2,  which  is  comparable  to  the  Kolmogorov  scale.  Because  of  the  higher 
density,  the  recombination  reactions  have  sped  up  and  are  in  between  the  two  turbulence  scales.  The 
thermal  NO*  and  CH2  +  N2  reactions  are  slow  relative  to  both  turbulence  scales,  but  CH  +  N2  is 
fast.  Time  scales  for  the  CH4  +  M  and  the  CO  +  OH  reactions  are  comparable  to  the  integral  and 
Kolmogorov  scales,  respectively. 

It  is  clear  that  while  some  combustion  reactions  fall  in  the  flamelet  regime  and  some  in  the  distri¬ 
buted  regime,  many  important  reactions  fall  in  between.  Turbulence-chemistry  interactions  occur  in  a 
greater  variety  of  modes  (Damkohler  numbers)  than  accounted  for  by  modern  combustion  theory. 
Most  modern  theories  classify  a  flame  as  lying  entirely  in  some  particular  domain  on  a  set  of  intensity 
or  time-scale  axes,  missing  the  critical  issue  that  there  are  multiple  interaction  scales  present  in  a  sin¬ 
gle  flame. 

The  simultaneous  prediction  of  kinetically  sensitive  phenomena  such  as  pollutants,  flame  stability 
and  combustion  efficiency  will  require  recognition  of  the  multiplicity  of  interactions  within  the  same 
flame. 
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Table  I 

Characteristic  Time  Scales  of  Selected  Reactions  in  Two  Combustors 


Reactions 

0.1  atm 

H2  scramjet 
(Tc,  s) 

10  atm 

CH4  combustor 
(Tc,  s) 

1. 

H  +  02  =  OH  +  0 

3.49  xlO'5 

1.67  xlO'7 

2. 

O  +  H2  =  OH  +  H 

2.64  xlO'5 

1.81  x  10'5 

3. 

H  +  OH  +  M  =  H20  +  M 

1.59  xlO'3 

7.20  xlO-5 

4. 

N2  +  0  =  NO  +  N 

9.45  xlO4 

3.92  x  10  ~2 

5. 

CH4  +  M  =  CH3  +  H  +M 

4.68  xlO"4 

6. 

CH4  +  OH  =  CH3  +  H20 

2.48  x  10'7 

7. 

ch3  +  OH  =  ch2  +  h2o 

1.85  xlO'6 

8. 

CO  +  OH  =  C02  +  H 

2.22  xlO'5 

9. 

CH  +  N2  =  HCN  +  N 

3.43  xlO"6 

10. 

CH2  +  N2  =  HCN  +  NH 

0.20 

Table  II 

Characteristic  Turbulence  Time  Scales 


u'  (m/s) 

At  (m) 

»■/(  s) 

e  (m2/s3) 

r*(  s) 

Scramjet  combustor 

300 

0.002 

6.67  xlO'6 

1.35  xlO10 

6.01  xlO'7 

Turbine  combustor 

10 

0.005 

5.0  xlO'4 

2.0  x  10s 

1.81  x  10'5 
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Section  3 
CONCLUSIONS 

The  goal  of  this  three-year  research  program  was  a  quantitative  understanding  of  turbulence- 
chemistry  interactions  in  turbulent  gaseous  flames  under  conditions  approaching  blowoff.  The  pro¬ 
gram  built  on  prior  work,  first  on  unassisted  and  more  recently  on  pilot-stabilized  CO/H2/N2  jet 
flames  in  coflowing  air. 

Although  the  fast  chemistry  or  “mixed-is-burned”  paradigm  is  useful  for  understanding  traditional 
design  issues,  present-day  demands  on  combustion  equipment  and  thus  on  computational  models  are 
increasing  the  need  for  a  more  advanced  understanding  of  turbulence-chemistry  interactions.  For  ex¬ 
ample,  (i)  flameout  and  relight  in  turbine  combustors  are  related  to  interactions  of  turbulence  with 
chain-branching  reactions,  (ii)  hydrogen  burnout  in  supersonic  combustors  is  related  to  interactions 
with  recombination  reactions,  and  (iii)  emissions  of  NO*,  CO,  smoke  and  other  observables  are  relat¬ 
ed  to  nonequilibrium  among  species  such  as  oxyhydrogen  radicals  and  Q.Hy. 

The  main  technical  conclusions  of  this  work  are 

1.  To  obtain  higher  turbulent  strain  rates  than  attainable  with  jet  flames,  a  bluff-body  stabilized 
burner  was  used.  It  will  be  necessary  to  accurately  model  such  flows,  which  are  more  complex 
as  regards  both  turbulence  and  numerics  than  jet-like  flows,  if  the  link  between  theory  and  ex¬ 
periment  is  to  continue  as  fruitfully  as  before. 

2.  Although  originally  developed  for  H2  flames,  Raman  systems  can  be  modified  to  operate  use¬ 
fully  even  in  the  sooting/chemiluminescent  environment  of  hydrocarbon  flames.  Measurements 
were  made  in  CO/H2/N2  and  CH,  flames. 

3.  Complementarity,  the  bluff-body  stabilized  burner  was  modeled: 

i.  A  thermochemical  submodel  based  on  partial  equilibrium  in  the  oxyhydrogen  radical  pool 
was  developed  for  the  27.5%  CO/  32.3%  H 2/  40.2%  N2-air  system.  The  chemistry  can 
be  described  in  terms  of  two  scalars. 

ii.  The  elliptic  form  of  the  time-averaged  Navier  Stokes  equations  with  k-e  closure  was 
solved  using  an  iterative  finite  volume/pressure  correction  algorithm.  Mean  properties 
such  as  density  were  obtained  at  each  iteration  by  convolution  with  the  joint  pdf  over  the 
two  thermochemical  scalars.  The  pdf  itself  was  obtained  by  the  moment 
equation/assumed  shape  method. 

4.  Calculations  were  compared  with  the  data.  Agreement  was  very  encouraging  on  mixture  frac¬ 
tion  mean  and  variance,  temperature,  and  species  concentration  fields.  Numerical  error  in  the 
mean  hydrodynamic  field  calculations  can  be  demonstrably  reduced  to  acceptable  levels.  The 
k-e  turbulence  model  gave  too  rapid  an  initial  decay.  To  improve  the  turbulence  model  and  to 
have  a  formalism  permitting  three  or  more  scalars  as  required  for  hydrocarbon  fuels,  the  pdf 
transport  method  should  be  merged  with  conventional  solvers  for  the  mean  hydrodynamics. 

5.  Preliminary  work  on  Monte-Carlo  pdf  transport  methods  in  recirculation  stabilized  flames. 

Second-order  closure  of  the  turbulence  correlations  could  have  been  employed  in  this  study 
to  relax  some  of  the  k-e  model  assumptions.  The  plan,  however,  is  to  proceed  directly  to  the 
Monte-Carlo  joint  pdf  transport  method,  which  is  closed  to  all  orders  and  will  provide  the  same 
level  of  modeling  sophistication  in  recirculating  flow  as  already  demonstrated  for  jet  flames 
[Pope  and  Correa  (1987)]. 
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6.  The  challenge  of  developing  universal  turbulence-chemistry  models  can  be  assessed  by  direct 
calculation  of  a  matrix  of  Damkohler  numbers,  Dy=  r,/r;,  where  j  represents  the  chemical  reac¬ 
tion  of  interest  and  i  the  turbulence  scale  of  interest.  It  becomes  apparent  that  while  some 
combustion  reactions  fall  in  the  flamelet  regime  and  some  in  the  distributed  regime,  many  im¬ 
portant  reactions  fall  in  between.  Turbulence-chemistry  interactions  occur  in  a  greater  variety 
of  modes  than  accounted  for  by  modem  combustion  theory.  The  simultaneous  prediction  of 
kinetically  sensitive  phenomena  such  as  pollutants,  flame  stability,  and  combustion  efficiency 
will  require  a  new  approach. 

Suggestions  for  future  work  are 

i.  Laser  Doppler  Velocimetry,  minor  species  measurements  (laser-induced  fluorescence),  and 
joint  major-minor  species  measurements  (LIF-spontaneous  Raman)  in  recirculation-stabilized 
flames. 

ii.  Development  and  assessment  of  the  hybrid  Navier-Stokes-for-mean-flow/pdf-  transport 
method,  as  outlined  above. 
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Section  5 

WRITTEN  PUBLICATIONS 

The  following  publications  were  made  during  the  reporting  period  May  1988  -  May  1991,  support¬ 
ed  wholly  or  in  part  by  this  contract: 

1.  Gulati,  A.  and  Correa,  S.M.,  “Local  Extinction  Due  to  Turbulence  in  Non-premixed  Flames,” 
AIAA  Paper  No.  87-1717, 23rd  AIAA  Joint  Propulsion  Conference,  San  Diego,  CA,  June  29- 
July2, 1987. 

2.  Correa,  S.M.  and  Gulati,  A.,  “Non-premixed  Turbulent  CO/H2  Flames  at  Local  Extinction 
Conditions,”  Twenty-Second  (International)  Symposium  on  Combustion,  The  Combustion  In¬ 
stitute,  Pittsburgh,  PA,  pp.  599-606, 1988. 

3.  Gulati,  A.  and  Correa,  S.M.,  “Raman/LV  Measurements  and  Modeling  in  a  CO/H2/N  2  Flame 
at  High  Reynolds  Number,”  Paper  89-0056, 27th  AIAA  Aerospace  Sciences  Meeting,  Reno, 

NV,  Jan.  9-12, 1989. 

4.  Correa,  S.M.,  “Appropriate  Turbulence-Chemistry  Models  for  Nonpremixed  Combustion,” 
Chemical  and  Physical  Processes  in  Combustion,  Class  I  Report  89CRD179,  GE  Corporate 
Research  and  Development  Center,  Schenectady,  New  York,  1989,  also  Eastern  States  Section 
of  Combustion  Institute,  Fall  Technical  Meeting,  Albany,  NY,  Oct.  30  -  Nov.  1, 1989. 

5.  Correa,  S.M.,  “Relevance  of  Non-premixed  Laminar  Flames  to  Turbulent  Combustion,” 
Proceedings  of  NASA  Langley/ICASE  Combustion  Workshop,  Oct.  2-4, 1989,  to  be  published 
by  Springer-Verlag. 

6.  Correa,  S.M.,  “A  Review  of  NO*  Formation  Under  Gas-Turbine  Combustion  Conditions,”  sub¬ 
mitted  to  Combustion  Science  and  Technology,  December  1990. 

7.  Correa,  S.M.  and  Gulati,  A.,  “Modeling  Measurements  in  a  Bluff-Body  Stabilized  Flame,”  sub¬ 
mitted  to  Combustion  and  Flame,  August  1991. 
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Section  6 

PROFESSIONAL  PERSONNEL 

Dr.  Sanjay  M.  Correa  and  Dr.  Anil  Gulati  were  principal  investigators  on  this  project.  They  were 
supported  by  Mr.  Frank  Haller,  laser  diagnostics  specialist,  and  Ms.  Janet  Sober,  computer  analyst. 

Dr.  Correa  earned  his  Ph.D.  (Aerospace  Engineering)  in  1981  from  the  University  of  Michigan, 
with  a  thesis  entitled,  “Group  Behavior  of  Liquid  Fuel  Sprays.”  Dr.  Gulati  earned  his  Ph.D. 
(Aerospace  Engineering)  in  1986  from  the  University  of  Michigan,  with  a  thesis  entitled,  “Simultane¬ 
ous  Density-Velocity  Measurements  in  Premixed  Turbulent  Flames  to  Assess  a  Theoretical  Model.” 

Section  7 
INTERACTIONS 

The  following  external  interactions  occurred  or  were  arranged  during  the  reporting  period  May 
1988  to  June  1991: 

University  Seminars 

Ohio  State  University  (January  1989) 

University  of  Connecticut  (November  1989) 

Cornell  University  (December  1989) 

MIT  (November  1990) 

Invited  Talks 

Current  Problems  in  Gas-Turbine  Combustion 
Fall  Technical  Meeting 

Eastern  States  Section  of  the  Combustion  Institute 
December  3-5, 1990 
Orlando,  Florida. 

Successes  and  Failures  of  Numerical  Combustion 

SIAM  Fourth  International  Symposium  on  Numerical  Combustion 

December  2-4, 1991 

St.  Petersburg,  Florida. 

Other  Interactions  in  which  AFOSR-Sponsored  Research  was  Discussed 

SANDIA  Turbulent  Combustion  Research  Group  meetings  (2/year) 

AIAA  Short  Course,  Leader:  Computational  Methods  in  Combustion:  Gas-Turbines,  Ramjets 
and  Scramjets,  Orlando,  FL,  July  1990. 
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